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ABSTRACT 


Lamberson,  Steven  Edward.  Ph.D.,  Purdue  University,  August  1985. 
Equivalent  Continuum  Finite  Element  Modelling  of  Plate-like  Space 
Lattice  Structures.  Major  Professor:  Henry  Yang. 

A  variety  of  research  projects  are  being  pursued  involving  the 
dynamics  and  control  of  large  plate-like  sf^ce  platforms  made  up 
of  repetitive  lattice- type  truss  structures.  A  method  involving 
finite  element  modelling  of  an  equivalent  continuum  formulation 
based  on  matching  the  strain  energy  and  kinetic  energy  is  developed 
for  truss  type  lattices  with  pinned  joints.  The  method  is  shown  to 
give  modal  results  consistent  with  those  obtained  using  detailed 
finite  element  modelling  of  the  pin  jointed  space  lattice  structure, 
even  for  structures  with  fairly  small  numbers  of  repetitions  of  an 
identical  unit  cell.  Feedback  controllers  designed  using  reduced 
system  models  derived  from  these  modal  results  using  Modal  Cost 
Analysis  are  shown  to  perform  as  well  as  controllers  designed  using 

the  detailed  analysis  results.  _ 

/ 

The  efficiency  of  this  method  for  coupled  structure,  control 
system  design  is  demonstrated  with  a  parametric  study.  Lattice  plate 
finite  elements  are  used  to  examine  the  effect  of  variation  of 
several  fundamental  structural  parameters  on  the  natural  frequencies 
and  mode  shapes  of  the  structure.  Feedback  control  systems  are 

i 

designed  and  the  resulting  system  performance  evaluated  by 


examining  the  steady  state  regulation  cost  of  the  structure  as  a 
function  of  the  structural  design  parameters. 


A  micropolar  plate  continuum  model  of  large  plate-like  repetitive 
space  lattice  structures  with  rigid  joints  is  derived.  A  plate  finite 
element  is  derived  based  on  this  continuum  model  with  micropolar 
rotations  and  transverse  shear  deformations  included  as  nodal  degrees 
of  freedom  ~( d . o . f .  4~s')>.'  The  natural  frequencies  and  mode  shapes  are 
calculated  using  this  element  for  a  free  floating  hexahedral  plate-like 
space  lattice  structure.  These  natural  frequencies  and  mode  shapes 
are  compared  to  those  calculated  using  a  detailed  finite  element 
model  (with  every  structural  member  modelled  by  a  beam-column 
element)  for  several  sets  of  frame  member  cross  sectional  properties. 
The  static  deflections  of  a  centrally  loaded  and  corner  supported 
rectangular  plate-like  space  lattice  structure  are  also  obtained  and 
compared  using  these  two  types  of  structural  models. 
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CHAPTER  I  -  INTRODUCTION 


Problem  Statement 

Now  that  the  shuttle  has  made  access  to  space  more  economical, 
many  projects  are  being  proposed  which  involve  large  space  platforms 
made  up  of  lattice- type  truss  structures.  These  structures  are 
generally  too  flexible  to  meet  the  mission  requirements  and  automated 
control  systems  must  be  used  both  to  maintain  the  structures'  orienta 
tion  and  to  reduce  vibrations  of  the  structure.  Often,  these  large 
space  structures  are  made  up  of  a  simple  cell  of  truss  bars  repeated 
a  large  number  of  times  in  one  or  two  directions.  The  resulting 
structure  resembles  a  beam  or  a  plate  when  viewed  from  a  distance. 
Detailed  finite  element  models  of  such  truss  structures  have  a  large 
number  of  degrees  of  freedom  (d.o.f.'s),  which  are  cumbersome  and 
expensive  to  use  in  modal  analysis  even  when  dynamic  substructuring 
and  reduction  methods  are  used  extensively.  An  alternative  method 
is  to  derive  a  plate  finite  element  model  based  on  an  equivalent 
continuum  formulation  which  is  used  to  perform  the  modal  analysis. 

In  either  case  the  modal  model  generated  must  be  reduced  to  an 
evaluation  model  of  a  size  which  can  be  used  to  evaluate  the 
effectiveness  of  various  control  systems.  The  modal  model  must  also 
be  reduced  further  to  a  series  of  design  models  which  are  used  to 
design  feedback  control  systems.  For  large  lattice-type  space  truss 


structures  the  order  of  the  detailed  finite  element  model,  the  order 
of  the  solvable  eigenvalue  problem,  the  order  of  the  solvable  control 
system  evaluation  model,  and  the  order  of  the  controller  which  can 
be  implemented  usually  from  a  rapidly  decreasing  monotonic  sequence. 
The  object  of  this  research  is  to  develop  equivalent  plate  finite 
elements  based  on  equivalent  continue  representations  of  large 
plate-like  space  lattices.  The  resulting  equivalent  plate  finite 
elements  are  shown  to  be  accurate,  flexible,  and  efficient.  They 
are  also  demonstrated  to  be  useful  in  designing  and  evaluating  the 
closed  loop  system  performance  of  large  repetitive  plate-like  space 
lattices  with  reduced  order  controllers. 

Structural  Modelling  Methods 
Background 

There  have  been  many  studies  examining  finite  element  modelling 
methods  of  various  types  of  structures  which  will  yield  accurate 
modal  results.  There  have  also  been  a  variety  of  methods  developed 
to  allow  progressively  larger  finite  element  models  to  be  used  for 
modal  analysis.  In  general,  this  is  done  using  some  means  of  model 
reduction.  Several  methods  have  been  developed  to  allow  large  finite 
element  models  (which  are  discrete  representations)  to  be  reduced  to 
smaller  mathematical  models  (which  are  also  discrete  representations) 
These  methods  range  from  heuristic  techniques  such  as  Guyan 

reduction*  which  rely  heavily  on  the  user's  experience  for  accuracy, 

2 

to  iterative  techniques  such  as  generalized  dynamic  reduction  ,  which 


effectively  eliminate  the  need  for  an  experienced  user  at  the  expense 


of  substantially  increasing  the  computational  cost.  Other  techniques 
involve  breaking  the  structure  up  into  substructures  which  are 
analyzed  dynamically  using  a  system  modelling  technique  such  as 
Component  Modal  Synthesis.  Basically,  this  involves  developing 
a  modal  model  of  each  of  the  subcomponents,  then  combining  these 
individual  modal  models  into  a  system  model  whose  only  physical 
d.o.f.'s  are  the  boundaries  connecting  the  various  substructures. 

The  remaining  d.o.f.'s  are  those  representing  the  modal  coordinates 
of  the  individual  substructures.  The  principle  problem  with  this 
method  is  that  it  is  often  difficult  to  determine  which  modes  of  a 
substructure  will  contribute  to  the  modes  of  the  entire  structure 
which  are  within  the  desired  frequency  bandwidth.  Therefore,  while 
the  available  techniques  have  increased  the  order  of  finite  element 
models  of  structures  that  can  be  utilized  to  perform  structural 
vibration  analysis,  there  are  sufficient  difficulties  associated 
with  each  reduction  method  to  require  that  the  original  finite 
element  model  of  the  system  be  generated  using  the  minimum  number  of 
d.o.f.'s  consistent  with  generating  an  accurate  solution4. 

Since  many  of  the  large  flexible  space  structures  being 
considered  are  built  on  a  platform  which  is  a  space  lattice 
consisting  of  truss  type  members  connected  together  in  a  repetitive 
manner,  several  methods  have  been  developed  to  allow  the  repetitive 
space  lattice  to  be  modelled  as  an  equivalent  continuum.  The  first 
of  these  applies  static  loads  to  a  model  of  a  unit  cell  of  the 
repetitive  space  lattice  in  such  a  way  that  the  displacement 
calculated  can  be  used  to  determine  one  or  more  of  the  "equivalent" 
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properties.  Another  method  is  to  calculate  the  equivalent  properties 
based  on  the  number  of  truss  elements  within  the  unit  cell  aligned 
in  each  direction.6  A  third  method  matches  the  strain  energy  of  the 
unit  cell  in  terms  of  a  Taylor's  series  expansion  of  the  midplane 
displacements  with  the  strain  energy  of  the  unit  cell  in  terms  of  a 
Taylor's  series  based  on  the  midplane  strains.^  The  existing  work 
has  focused  on  predicting  the  vibration  modes  of  simply  supported 
truss  lattices  with  pinned  joints.  There  has  been  some  work  with 
simply  supported  "beam-like"  lattices  with  the  frame  members  connected 
by  rigid  joints  which  makes  use  of  the  higher  order  micropolar  beam 
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continuum  formulation. 


Detailed  Models 

For  this  research  the  baseline  against  which  the  performance  of 
the  equivalent  continuum  models  is  evaluated  are  calculated  values 
of  the  free  vibration  natural  frequencies  and  mode  shapes.  For 
repetitive  plate-like  space  lattices  with  pinned  joints  the  detailed 
finite  element  model  is  made  with  each  truss  member  modelled  using 
a  single  axial  force  type  truss  element.  Chapter  II  contains  a 
discussion  of  the  detailed  finite  element  modelling  method,  a 
description  of  the  detailed  finite  element  model  of  a  specific  space 
lattice,  and  a  comparison  of  the  resulting  natural  frequencies  and 
mode  shapes  with  those  obtained  with  an  equivalent  transverse  shear 
plate  finite  element.  Repetitive  plate-like  space  lattices  with 
rigid  joints  are  modelled  using  a  detailed  finite  element  model  with 
each  frame  member  modelled  as  a  beam-columm  type  finite  element. 


Chapter  IV  contains  a  description  of  the  detailed  finite  element 
model  of  a  specific  space  lattice,  and  a  comparison  of  the  resulting 
natural  frequencies  and  mode  shapes  with  those  obtained  with  an 
equivalent  micropolar  plate  finite  element.  In  both  cases,  Guyan 
reduction  is  used  to  reduce  the  order  of  the  eigenvalue  problem. 

Equivalent  Modelling  of  Space  Lattices  with  Pinned  Joints 
The  equivalent  plate  finite  element  for  repetitive  space  lattices 
with  pinned  joints  is  derived  using  an  existing  equivalent  continuum 
representation. ^  The  strain  energy  expression  derived  in  the 
reference  in  terms  of  the  in-plane,  bending,  and  transverse  shear 
strains  is  used  as  a  constitutive  law  to  develop  a  transverse  shear 
type  plate  finite  element.  This  element  assumes  a  linear  displacement 
functions  for  the  two  inplane  displacement  and  two  transverse  shear 
deformation  d.o.f.'s  at  each  node,  and  a  bicubic  Hermite  polynomial 
displacement  function  for  the  transverse  displacement.  This  element 
allows  considerable  simplicity  and  flexibility  in  handling  various 
geometries,  boundary  conditions,  attachments,  and  other  structural 
complexities.  Chapter  II  contains  a  discussion  of  the  existing 
equivalent  continuum  representation,  a  derivation  of  the  equivalent 
plate  finite  element,  a  description  of  the  equivalent  plate  finite 
element  model  for  a  specific  space  lattice,  and  a  comparison  of  the 
resulting  natural  frequencies  and  mode  shapes  with  those  obtained 
using  a  detailed  finite  element  model.  Chapter  III  contains  a 
parametric  study  examining  the  effect  on  the  natural  frequencies 
and  mode  shapes  of  varying  the  geometric  properties  of  the  unit  cell 
building  block  of  a  specific  repetitive  space  lattice. 


Equivalent  Modelling  of  Space  Lattices  with  Rigid  Joints 

An  equivalent  continuum  representation  is  available  for  beam-like 
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repetitive  space  lattices  with  rigid  joints.  This  method  is 
extended  to  allow  plate-like  repetitive  space  lattices  with  rigid 
joints  to  be  represented  by  a  micropolar  plate  continuum.  The  strain 
energy  expression  derived  for  this  continuum  representation  in  terms 
of  the  in-plane,  bending,  transverse  shear,  and  micropolar  strains 
is  used  as  a  constitutive  law  to  develop  a  micropolar  type  plate 
finite  element.  This  element  assumes  linear  displacement  functions 
for  the  two  in-plane  displacement,  two  transverse  shear  deformation, 
and  two  micropolar  rotation  d.o.f.'s  at  each  node,  and  a  bicubic 
Hermite  polynomial  displacement  function  for  the  transverse  displace¬ 
ment.  Chapter  IV  contains  a  derivation  of  the  micropolar  equivalent 
continuum  and  micropolar  type  plate  finite  element  of  repetitive 
space  lattices  with  rigid  joints.  Chapter  IV  also  contains  a 
description  of  the  micropolar  plate  finite  element  model  for  a 
specific  space  lattice  and  a  comparison  of  the  resulting  natural 
frequencies  and  mode  shapes  for  several  sets  of  cross  sectional  areas 
with  those  obtained  using  detailed  finite  element  models. 

It  should  be  pointed  out  that  by  using  the  equivalent  continuum 
representations  as  bridges  between  the  discrete  space  lattice 
structures  and  the  discrete  equivalent  plate  finite  element  models, 
the  number  and  location  of  node  points  in  the  equivalent  plate  models 
is  determined  based  on  the  number  of  modes  which  are  desired  to  be 
accurately  calculated,  not  by  the  geometry  of  the  structure. 
Specifically,  there  is  no  requirement  that  (or  advantage  to)  having 
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the  internal  nodal  points  of  the  equivalent  plate  finite  element 
models  coincide  with  joint  locations  or  unit  cell  boundaries.  For 
large  structures,  the  number  of  nodal  d.o.f.'s  used  in  the  equivalent 
plate  finite  element  model  will  generally  be  much  smaller  than  the 
number  of  nodal  d.o.f.'s  required  to  represent  the  structure  as  a 
detailed  finite  element  model. 

Control  System  Design  Method 

Even  though  the  order  of  the  eigenvalue  problem  describing  the 

repetitive  plate-like  space  lattice  structure  has  been  reduced 

significantly  using  Guyan  reduction  for  the  detailed  finite  element 

models  and  the  inherent  reduction  associated  with  using  the  equivalent 

plate  finite  elements,  there  will  still,  in  general,  be  many  more 

natural  frequencies  and  mode  shapes  extracted  than  can  be  readily 

used  to  evaluate  control  system  performance.  In  addition,  the  higher 

natural  frequencies  and  mode  shapes,  particularly  when  Guyan  reduction 

g 

is  used  are  known  to  contain  considerable  error.  Therefore,  the 
higher  natural  frequencies  and  mode  shapes  are  truncated  and  those 
retained  are  used  as  a  modal  model  of  the  structure  for  evaluation 
purposes.  In  order  to  evaluate  the  performance  of  the  closed  loop 
system  containing  this  modal  model  of  the  structure  and  a  feedback 
controller  a  cost  function  is  defined  as  the  integral  over  all  time 
of  the  weighted  norm  of  a  vector  of  system  d.o.f.'s  of  interest.  The 
optimum  structure,  feedback  controller  combination  is  defined  to  be 
that  system  which  minimizes  this  quadratic  cost  function.*®  In 
general,  the  order  of  the  controller  to  be  used  is  much  smaller 
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than  the  order  of  the  modal  model  needed  to  represent  the  dynamic 
properties  of  the  structure.  Therefore,  the  modal  evaluation  model 
must  be  further  reduced  to  obtain  the  modal  model  for  which  the 
optimal  control  system  will  be  designed.  While  this  could  also  be 
done  by  truncating  the  higher  modes,  an  alternative  method  called 
Modal  Cost  Analysis  (MCA)  is  used  instead.**  In  this  method  the 
open  loop  cost  of  each  mode  is  used  to  rank  the  modes  in  terms  of 
their  effect  on  system  performance.  Those  modes  which  make  the 
largest  contribution  to  the  open  loop  system  cost  are  retained.  Once 
the  reduced  order  control  system  design  modal  model  has  been  obtained 
the  standard  Linear  Quadratic  Gaussian  (LQG)  theory*0  is  used  to 
design  a  feedback  controller  which  is  optimal  for  the  reduced  order 
system.  This  system  performance  is  then  calculated  using  the  reduced 
order  controller  to  drive  the  evaluation  modal  model.  It  is  important 
that  the  evaluation  modal  model  be  of  substantially  higher  order  than 
the  reduced  order  design  model  to  determine  the  effect  of  control 
system  spillover  on  the  modes  which  were  not  retained  in  the  reduce 
order  model.  Chapter  II  contains  a  description  of  the  feedback 
control  system  design  and  closed  loop  system  performance  evaluation 
methods  used.  This  chapter  also  contains  a  description  of  the 
application  of  these  methods  for  modal  models  of  a  specific  lattice 
structure  generated  with  both  a  detailed  finite  element  model  and  an 
equivalent  plate  finite  element  model,  and  a  comparison  of  the 
resulting  closed  loop  system  performance  costs  for  several  different 
orders  of  reduced  order  controller.  Chapter  III  describes  the 
application  of  these  control  system  design  and  evaluation  methods  to 


a  series  of  modal  models  generated  using  equivalent  plate  finite 
element  models  with  various  member  properties  in  the  form  of  a 
parametric  study.  The  resulting  system  performance  costs  are 
used  to  determine  the  best  set  of  member  properties. 


CHAPTER  II  -  EQUIVALENT  MODELLING  OF 
TRUSSES  WITH  PINNED  JOINTS 


Many  large  space  structures  are  being  designed  which  are  made 

up  of  a  simple  cell  of  truss  bars  repeated  a  large  number  of  times 

in  two  directions.  Such  a  plate-like  space  lattice  has  been 

represented  as  an  equivalent  plate  continuum  with  transverse  shear 

deformation  included. ^  This  representation  is  used  to  derive  an 
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equivalent  plate  finite  element.  As  a  benchmark  to  compare  the 
resulting  natural  frequencies  and  mode  shapes  of  the  structure  to, 
a  detailed  finite  element  model  is  also  developed  in  which  each 
member  of  the  truss  is  modelled  as  a  simple  axial  force  type  truss 
element.  A  specific  plate-like  space  lattice  structure  is  used  to 
demonstrate  the  accuracy  of  the  simplified  equivalent  plate  finite 
element  model.  The  eight  by  eight  repeating  cell  structure  chosen 
is  large  enough  to  allow  the  continuum  modelling  to  be  acceptable. 
However,  most  structures  would  have  many  more  repetitions  than  this 
example.  This  would,  of  course,  increase  the  cost  and  cumbersomeness 
of  the  detailed  finite  element  model  relative  to  the  continuum  model. 
This  would  also  improve  the  basic  assumptions  inherent  in  the  continuum 
model  development.  The  results  shown  that  even  with  as  few  as  eight 
repeating  cells  in  each  direction,  the  control  system  designs  based 
on  the  simplified  finite  element  continuum  models  are  effective  as 
reduced  order  controllers  for  the  detailed  finite  element  model. 


***,*."  v*  V" .  »*,  V,.  >•",  *•*’-  { J"v  ■  ** 

»  .1  kJi  lA if  fc.i  M fc  1*3  *-A. H..V  JL.TL ft.."  ^  ■**-» 
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Detailed  Finite  Element  Model 

A  large  finite  element  model  of  the  lattice-type  truss  structure 

13 

is  developed  using  standard  axial  force  truss  bar  elements.  This 
model  treats  each  member  of  the  structure  as  a  pinned-pinned  axial 
force  member  with  a  constant  strain.  In  order  to  calculate  the 
frequencies  and  modes  of  the  detailed  finite  element  model,  the 
model  must  be  reduced  to  a  manageable  order  using  some  form  of 
dynamic  reduction.  Guyan  reduction*  is  used  in  this  study  to  reduce 
the  detailed  finite  element  model  to  a  dynamic  analysis  model  retain¬ 
ing  d.o.f.'s  specified  by  the  user.  They  are  normally  chosen  to 
adequately  represent  the  mass  distribution  within  the  structure  and 
the  anticipated  frequencies  and  mode  shapes. 

[M]{q}  +  [K]{ q}  =  0  (II. 1) 


where: 

{ q }  =  vector  of  the  d.o.f.'s  of  the  complex  model; 

[M]  =  mass  matrix  of  the  complex  model;  and 
[K]  =  stiffness  matrix  of  the  complex  model. 

The  mass  and  stiffness  matrices  are  partitioned  into  submatrices 
associated  with  the  d.o.f.'s  to  be  retained  and  those  to  be  removed. 


{ q_ }  =  vector  of  the  d.o.f.'s  to  be  retained; 

a 

{qQ}  =  vector  of  the  d.o.f.'s  to  be  removed;  and 
{F}  =  force  vector. 

The  static  problem  can  be  reduced  exactly  by  using  the  lower  partition 
of  Equation  (II. 3)  as  a  constraint  equation. 

'  l[Kaa]  -  =  (F1  <n'4 

Guyan  reduction  assumes  that  this  same  transformation  can  be  applied 
to  the  mass  matrix. 

^  -  W.03C,!ooT1t«,.]T  - 

+  <u'5 

In  general  this  reduction  is  not  exact.  However,  it  is  widely  used 
and  gives  good  results  if  a  suitable  set  of  d.o.f.'s  is  chosen  to 
be  retained. 

In  order  to  perform  free  vibration  analysis  of  a  lattice  plate 
structure  with  all  edges  free,  the  rigid  body  motion  of  the  structure, 
which  renders  the  stiffness  matrix  singular,  must  be  accounted  for. 

In  this  study,  the  procedures  used  in  the  NASA  Structural  Analysis 
code  (NASTRAN)  will  be  used.  This  procedure  is  given  in  detail  in 
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the  NASTRAN  theoretical  manual  and  is  briefly  described  here.  The 
reduced  dynamic  system  is  solved  for  the  eigenvalues  and  eigenvectors. 

{[Kl  -  X [M  ] }  {q  (a)}  =  {0}  (II. 6) 

a  a  a 

where: 

a.  =  the  ith  eigenvalue  of  the  reduced  dynamic  system; 

(q  .}  =  {q  (A.)}  =  the  ith  eigenvector  of  the  reduced 

dynamic  system;  and 

[Q  ]  =  [qa}»  qa2»---]  =  the  reduced  dynamic  system  modal 

matrix. 

However,  for  the  free  floating  or  unsupported  plate  the  stiffness 
matrix  is  singular  to  order  six.  The  first  six  eigenvalues  are  not 
precisely  zero  due  to  numerical  error  in  the  eigenvalue  extraction 
process.  Consequently,  the  eigenvectors  calculated  are  not  precisely 
the  rigid  body  modes.  Therefore  a  procedure  is  implemented  to  replace 
these  six  eigenvalues  and  eigenvectors  with  six  rigid  body  modes. 

The  stiffness  and  mass  matrices  which  have  already  been  reduced 
using  Guyan  reduction  are  further  partitioned  based  on  a  set  of 
d.o.f.'s  specified  by  the  user  as  necessary  to  specify  the  rigid 
body  motions  of  the  structure  (qr>.  The  remaining  d.o.f.'s  are 
placed  in  a  set  {q^}. 


(II. 7) 


{q^}  =  the  d.o.f.  s  in  the  flexible  set;  and 

{ qr }  =  the  d.o.f.'s  used  to  define  the  rigid  body 
motion  of  the  structure. 

For  a  rigid  body  mode  the  equations  of  motion  can  be  reduced  to  a 
constraint  equation  relating  the  flexible  d.o.f.'s  to  the  rigid  body 
motion  d.o.f.'s. 

{[Ka]  -  u>2[Ma]}  { qfl }  =  {0}  (II 

For  rigid  body  modes,  u>  =  0,  thus  Eq.  (II. 8)  reduces  to 


[Kff]{qf>  +  [Kfr]{qr>  =  {0} 

(II 

{qf>  =  [0]{qr> 

(II 

where 

[D]  =  -  [Kff]"1[Kfr];  and 

w  =  natural  frequency  of  the  mode. 

By  matching  the  kinetic  energy  of  the  structure  in  the  rigid  body 
mode  and  using  the  constraint  Eq.  (1 1. 10)  to  eliminate  {q^}  the 
rigid  body  mass  matrix  is  defined. 


Substituting  Eq.  (II. 10)  for  {q^}  in  (II. 11)  and  multiplying  out 
Eq.  (II. 11)  gives 
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<$'  {qr}T[Mr]{qp}  =  {qr}T[D]T[Mff][D]{qr}  +  C qp> T[D]T[Mf r] { qr> 

+  {qr)T[Mfr]T[D]{qr}  +  { qr )T[Mrr] { qr> 

where  the  rigid  body  mass  matrix 

[Mr]  =  [D]T[Mff][D]  +  [D]T[Mfr]  +  [Mfr]T[D]  +  [Nr(J  .  (11.12) 

A  set  of  rigid  body  modes  is  generated  such  that  they  form  an  orthogo¬ 
nal,  normalized  set  with  respect  to  the  rigid  body  mass  matrix. 

[Qro]T[Mr][Qro]  =  [I]  (11.13) 

The  rigid  body  eigenvectors  are  calculated  for  the  flexible  d.o.f. 's 
and  augmented  to  the  rigid  body  modes  to  give  the  rigid  body  mode 
shapes  of  the  entire  set  of  d.o.f.'s. 

[D][Qr0] 

[Qro] 

The  original  modal  matrix  calculated  using  Eq.  (II. 6)  is 
partitioned  into  the  approximate  eigenvectors  and  accurate  flexible 
eigenvectors. 

[Qa]  -  [[Qar]  j  [Qaf]]  CQ]  -  [[Qao]  j  [Qaf]]  (0.15) 

The  rigid  body  partitioned  matrix  [Q  ]  is  discarded  and  replaced  by 
the  more  accurate  rigid  body  mode  shapes  calculated  using  Eq. 

(11.13).  The  resulting  modal  matrix  is  orthogonalized  and  normalized 
to  give  the  system  modal  matrix. 
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[Q]T[Ma][Q]  =  [I]  (11.16) 

Equivalent  Continuum  Representation 
The  method  developed  in  Reference  7  is  used  to  generate  strain 
energy  and  kinetic  energy  expressions  for  the  structure  in  terms  of 
strain  components  of  the  plate  at  the  midplane.  Briefly,  the 
displacements  are  assumed  to  vary  linearly  through  the  thickness  of 
the  plate  (Fig.  1). 

u(x,y,z)  =  u°(x,y)  +  z<t>x(x,y) 

v(x,y,z)  =  v°(x,y)  +  z<)>y(x,y)  (11.17) 

w(x,y,z)  =  w°(x,y)  +  ze|(x,y) 

where 

(x,y,z)  are  the  coordinates  of  a  point  within  the  plate; 

(u,v,w)  are  the  displacements  along  (x,y,z); 

(<t>  ,4>  )  are  the  rotations  about  the  (-y,x)  axes  at  the 
x  y 

midplane  (z  =  0); 

(u°,v°,w°)  are  the  displacements  along  (x,y,z)  at  the 
midplane;  and 

e3  =  cOz  at  the  midplane. 


The  axial  strain  in  the  truss  members  of  a  typical  unit  cell  is 
expressed  in  terms  of  strain  components. 


k 

e 


l  I  Hjk  s-k 

i=l  j=l  1  J 


(11.18) 
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where: 

e  =  axial  strain  in  the  k'th  truss  bar; 

k  k 

e •  •  =  strain  component  of  e  ;  and 

*  J 

i< 

i.  =  direction  cosines  of  the  member. 

These  are  then  expanded  using  a  Taylor's  series  in  terms  of  the 
strain  components  and  their  derivatives  about  some  arbitrary  origin 
within  the  unit  cell.  The  fact  that  the  forces  associated  with 
certain  of  the  strain  components  and  derivatives  are  zeros  is  used 
to  reduce  the  strain  energy  expression  to  an  expression  in  terms  of 
eight  strain  components  at  the  plate  midplane.  The  truss  cell  geometry 
used  in  this  study  is  one  of  those  analyzed  in  the  reference. 


where 


“cell  ‘Kell  IV(C*} 


Ucell  =  strain  energy  of  a  unit  cell; 
acell  =  cross  sectional  area  of  a  unit  cell; 

[W  ]  =  equivalent  constitutive  matrix; 

,  ,  _  r  0  0  0  0  0  0  0  0  J. 

e22*  e12’  Kll’  k22’  k12’  e13*  e23^  ’ 

0  _  1  ,au  .  _|v x . 

Ei  j  “  2  'ay  ax;’ 


(11.19) 


-  1  l  4.  A  } 

Ei3  '  2  (ax.  *xiJ  • 


Jr-  .4 
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Equivalent  Finite  Element  Model 

15 

The  16  d.o.f.  rectangular  plate  element  is  modified  to  include 
transverse  shear  terms  and  in-plane  terms  yielding  a  32  d.o.f. 
element  (Fig.  2).  The  interpolation  functions  are  assumed  as  linear 
except  for  the  out-of-plane  displacement  which  is  assumed  as  a  bi-cubic 
Hermite  polynomial. 

u  =  a^  +  a2x  +  agy  +  a^xy 

v  =  a5  +  a6x  +  a7y  +  a8xy 

2  2  3  2 

w  =  a9  +  a10x  +  any  +  a12*  +  a13xy  +  a14y  +  a15x  +  a16x  y 

+  a1?xy2  +  algy3  +  aigx3y  +  a2Qx2y2  +  a21xy3  (11.20) 

+  a  +  a  x^v3  +  a  x3v3 

a22x  y  a23x  y  a24x  y 
e13  =  a25  +  a26X  +  a27y  +  a28Xy 

c23  =  a29  +  a30x  +  a31y  +  a32Xy 

Substituting  Eq.  (11.20)  into  Eq.  (11.19)  gives  the  strain  energy 
expression  for  the  element. 


2acell 


Jarea^ 


{a}T[B]T[Weq][B]{a}dydx 


where 


{ea>  =  [B]  {a};  {a}  = 


(II.2D 


(II.21) 
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[B] 


B !  0  0  0 

0  E $2  B^  B^  » 

0  0  0  b5 


[Bj] 


[b2] 


[b3] 


[b4] 


[b5] 


“0  10  y  0  0 

0  0 

0  0  0  0  0  0 

1  x 

* 

_0  0  0.5  0.5x  0  0.5 

0  0 . 5x_ 

“0  0  0  -2  0  0  -6x 

-2y  0 

0  -6xy“ 

0  0  0  0  0  -2  0 

0  -2x  -6y  0 

_0  0  0  0  -1  0  0 

-2y  -2y  0  -3x2_ 

-2y2  0  -6xy2  -2y3 

-6xy3 

-2x2  -6xy  -2x3  -6x2y 

-6x3y 

> 

_-4x  -3y2  -6x2y  -6xy2 

-9x2y2_ 

“Oil  y  0  0 

0  0“ 

0  0  0  0  0  0 

1  x 

;  and 

_0  0  0.5  0. 5x  0  0.5 

0  0.5y_ 

“lxyxyOOO  0“ 

_0  0  0  0  1  x  y  xy_ 

Relating  the  32  displacement  function  coefficients  to  the  32  d.o.f.'s 
of  the  element  the  strain  energy  of  the  element  is  transformed  in 
terms  of  the  element  d.o.f.'s. 
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tz 
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el  2a 


1  tg)T 


cell 


AREA 


'  [T]T[B]T[Weq][B][T]dydx{g}  (11.22) 


where: 


{a}  =  [T](g}  ; 

{g}  =  {g2>T.  {93)T»  {g4>T!IT ;  and  AREA  is  the  area 

of  the  rectangular  plate  finite  element. 


This  defines  the  stiffness  matrix  of  the  element. 


[Ke,]  * 


’cell  JAREA 


[T]T[B]T[Weq][B][T]dydx 


(11.23) 


The  kinetic  energy  expression  derived  in  Ref.  7  is  used  to  generate  the 
lumped  mass  matrix  of  the  element. 


K'E-  =  2  acel 1 


m 


0 


0 

0 

m. 


mr 


0 

0 


m 


0  mr 


0 

m 


1 


1  w 
0  IDj 

0  0 
m2  0 
0  m0 


(11.24) 


where  nig,  m^,  and  are  equivalent  lumped  masses  of  the  space  truss 
defined  in  the  reference. 


[Mei]  - 


'M1 

0 

0 

0 


0  M. 


0  M, 


(11.25) 
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“  SYMMETRY 

0  0  mQ 

mj  0  0  m2 

0  0  0  m2 

0  0  0  0  0  0 

0  0  0  0  0  0  0 


Control  System  Design  Method 

The  detailed  finite  element  model  is  reduced  to  an  evaluation 

modal  model  which  can  be  readily  analyzed  to  determine  the  effectiveness 

of  the  various  control  designs.  The  modes  retained  in  the  evaluation 

model  are  arbitrarily  selected  to  be  the  lowest  frequency  modes  of 

the  complex  model  that  are  both  controllable  and  observable.  The 
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evaluation  model  is  used  in  the  real  modal  form. 


Wj  0  0. . .0 

0  W2  0...0 

0  0  W3. .  .0 


00  WnKQn  Bn 


{y}  -  [Cr  C2,  ...  Cn]  ]Q2  (  {z}  -  [Mj,  M2,  ...  Mn]  Jq2  ( +{v} 


(11.26) 
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where 


{u}  is  the  input  force  vector  to  the  system; 

(w)  is  the  input  disturbance  vector  to  the  system; 

{v}  is  the  measurement  disturbance  vector  of  the  system; 
{z>  is  the  measurement  vector  of  the  system; 

(y }  is  the  output  vector  of  the  system; 


(Q • }  *  with  q..  =  the  modal  coordinate  associated  with 


[BJ  = 


mode  i 


jiwith  {bn-  }T  =  the  modal  coefficients  of  the 


1  J  applied  force  for  mode  i; 


^  > with  {d.}  =  the  modal  coefficient  of  the 
1  J  applied  disturbance  vector  for  mode  i; 


=  [Ci  0]  with  {c. }  =  the  modal  coefficients  of  the 
output  vector  for  mode  i ; 

=  [m.j  0]  with  {m.. }  =  the  modal  coefficients  of  the 
measurement  vector  for  mode  i ; 


[wi]  = 


with  uk  =  the  natural  frequency  of  the 
mode  i ;  and 

C  =  the  modal  damping  coefficient. 


The  full  models  available  from  either  the  complex  model  or  the 
equivalent  plate  finite  element  models  usually  generate  too  many 
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.>  modes  to  be  used  directly  as  control  design  models.  Therefore,  some 

means  of  model  reduction  is  required.  Several  dynamic  reduction 
schemes  have  been  developed  to  al'ow  the  control  problem  to  influence 
the  model  reduction;  of  these  Modal  Cost  Analysis  (MCA)  is  the 
most  straightforward  to  use.^  In  order  to  use  MCA,  the  relative 
importance  of  a  set  of  outputs,  each  of  which  is  a  linear  combination 
of  d.o.f.'s  of  the  structure  is  used  to  build  a  positive  definite 
output  weighting  matrix  [Q].  The  cost  which  is  to  be  minimized  by 
the  control  system  is  then  defined  as  the  sum  of  a  weighted  norm  of 
the  output  vector  plus  a  weighted  norm  of  the  measurement  vector. 

V  =  Vy  +  PVz  (11.27) 

where 

Km  fT  t 

V  «  lim  E  {y}'[Q]{y)dt  ; 

y  J0 

rT  t  -i 

Vz  =  lim  E  {z}  [Z]  A{z}dt  ; 

X-XD  o 

E  =  expected  value;  and 
p  =  measurement  weight. 

The  weighting  matrix  used  for  measurements  is  the  inverse  of  the 
covariance  matrix  of  the  measurement  errors. 


[Z]  =  E[{ v}{v}T] 


(11.28) 
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For  small  values  of  modal  damping  this  can  be  reduced  to  an  uncoupled 
set  of  equations  as  derived  in  Ref.  17. 

Vi  =  — ^~3  ({ci)T[Q]{ci}  +  {mi}T[Z]"1{mi})({bi}T{bi})  (11.29) 

4?wi 

Control  design  models  of  order  n  are  derived  from  the  detailed 
finite  element  and  the  equivalent  plate  model  by  retaining  the  n/2 
modes  having  the  largest  modal  costs.  Since  rigid  body  modes  have 
infinite  modal  cost,  they  must  always  be  retained  in  the  reduced  models 
if  they  are  controllable  and  observable. 

Each  control  design  model  is  used  to  design  a  standard  LQG 
controller.^ 

(q)  =  [A] { q }  +  [B]{u}  +  [D]{w} 

(y)  =  [C]{q } 

(z)  =  [M] { q }  +  v  (11.30) 


where 


E{w}  =  E{ v }  =  0; 

E[ {w( x ) } {w( t ) }T]  =  [ W] 6 ( t— t ) ; 

E[{v(t) }{v(t) }T]  =  [V]6( t-x); 

E[{w(t) }{q(0) }T]  =  E[{v(t) }{q(0) }T]  =  0;  and 
E[{v(t) }{w(t) }T]  =  0  . 


Starting  with  a  reduced  system  in  the  form  of  Eq.  (11.26)  a  state 
estimator  is  designed  to  allow  a  feedback  controller  to  be  designed. 


{q }  =  [A] {q >  +  [B]{u}  +  [F][{z]  -  [M] {q }]  (11.31) 

Where  the  filter  matrix  is  defined  as  a  Riccati  equation. 

[f]  -  [PHxiVr1 

[0]  =  [P][A]T  +  [A][P]  -  [P][M]T[V]'1[M][P]  +  [0][W][D]T 

(11.32) 

Once  the  estimate  of  the  state  is  available  a  state  feedback  controller 
is  designed. 

(u)  =  [ G] { q }  (11.33) 

Where  the  feedback  gain  matrix  [G]  is  defined  as  a  Riccati  equation. 

[G]  -  -  1  [R]'1[3]T[K] 

[0]  .  Ck][a]  +  [A]T[K]  -  i  CK3[B]tR3'1CB]TCK3  +  [C]T[Q][C] 

(11.34) 


Closed  Loop  System  Performance 

The  full  order  controller  based  on  the  evaluation  model  can  be 
evaluated  directly,  as  in  Ref.  10. 


V 


=  V  +  pV 
y  M  u 


V 

y 


1  im  - 

T-KX,  1 


fT 


{y}T[Q]{y}dt  =  tr[C]T[Q][C]  [[X]  +  [P]] 


(11.35) 


o 


V  =  1  im  4"  {u}  [R]{u}dt  =  tr[G]  [R][G][X] 

u  T-*»  ■'o 

where 

CO]  =  [X]{[A]  +  [B][G]}T  +  {[A]  +  [B]  [G] }  [X  ]  +  [F][V][F]T 

(11.36) 

In  order  to  evaluate  the  full  evaluation  model  driven  by  a  reduced 
order  controller,  the  evaluation  model  must  be  augmented  with  the 
controller  model  as  in  Ref.  16. 

For  the  evaluation  model, 

{qe>  =  [Ae]{qe>  +  [Be](u)  +  [De]{w} 

{yi  =  tce]{qe} 

iz}  =  [Me]{qe>  +  {v}  (11.37) 

For  the  controller  model, 

{q>  =  [Ac]{q}  +  [Fr]{z> 

fu}  =  [GR] {q}  (11.38) 

where 

C^q]  =  [ar]  +  [BR][GR]  -  [Fr][Mr] 

This  forms  an  augmented  system. 
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{qa>  =  CAa]{qa>  +  [Dq]{wa)  (11.39) 


w 

C  0 

iV 

W  0 

• 

=  [C  3(q  >  [WJ  = 

M 

0  gr_ 

1 q 

a  a  a 

_o  v_ 

from  which 

Vy  =  tr[Ce]T[Q][Ce][Xe];  and 
vu  =  trEGR]T[R][GR3[xR] 


where 


and 


[o]  -  [Xa][Aa]T  +  [Aa][Xa]  +  [Da][Ua][Da]T  . 


Results 

Detailed  Finite  Element  Model 

A  unit  cell  of  the  lattice-type  space  structure  to  be  analyzed  in 

this  study  is  shown  in  Fig.  3.  The  dimensions  of  the  unit  cell  are 

those  from  Ref.  7.  It  is  15  meters  square  by  7.5  meters  deep. 

2 

Structural  members  on  the  upper  surface  have  an  area  of  80  mm  .  The 

2 

bottom  surface  members  have  an  area  of  50  mm  .  Elements  connecting 

o 

the  two  surfaces  have  an  area  of  10  mm  .  The  material  has  a  Young's 
Modulus  of  71.7x10^  Nt/rn^  and  a  density  of  2768  kg/m^.  The  structure 
to  be  analyzed  consists  of  identical  cells  repeated  eight  times  in 


I 


-  UPPER  SURFACE  MEMBERS 

-  LOWER  SURFACE  MEMBERS 

-  -  VERTICAL  MEMBERS 

-  DIAGONAL  MEMBERS 


ORIGINAL 
VALUES 
A  ,=80  mm2 

A^Omm2 
Av=l  Omm2 
A^=  I  Omm2 
L=  h=75m 


Figure  3.  Plate-like  Space  Lattice  Unit  Cell 


each  of  two  directions  (Fig.  4).  The  resulting  structure  resembles  a 
plate.  The  detailed  finite  element  model  is  made  using  a  pinned  truss 
element  for  each  member.  This  model  is  analyzed  using  NASTRAN.  The 
full  model  has  1254  d.o.f.'s.  The  reduced  dynamic  analysis  model  has 
214  d.o.f.'s  as  shown  in  Fig.  5.  This  structure  was  chosen  to  be 
large  enough  to  allow  an  equivalent  continuum  model  to  be  accurate 
over  a  significant  frequency  range,  but  small  enough  to  be  readily 
analyzed. 


Equivalent  Plate  Models 

The  truss  properties  of  the  unit  cell  shown  in  Fig.  3  are  used  to 

generate  the  strain  energy  and  kinetic  energy  of  the  structure  in 

terms  of  the  strain  components  at  an  arbitrary  point  within  the  cell. 

The  matrices  [W  ]  and  [Man]  in  the  strain  energy  expression  Eq.  (11.21) 
eq  eq 

and  kinetic  energy  expression  Eq.  (11.24),  respectively,  are  tabulated 
below  and  are  in  agreement  with  those  presented  in  Ref.  7. 

Constitutive  Matrix 


(11.39 


Hia888SB8*Sl!S8iei!S38*i0e*(S8»$3 
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Mass  Density  Matrix 


JV- 

kg/m 


0.17794 

0.0 

0.0 

0.14176m 

0.0 

0.17794 

0.0 

0.0 

0.14176m 

0.17794 

0.0 

0.0 

SYMMETRY 

2.3698m2 

0.0 

2.3698m2 

(11.40) 


These  expressions  are  used  to  derive  an  equivalent  plate  finite 
element.  This  element  is  used  to  form  two  equivalent  plate  models, 
a  coarse  one  (4x4)  containing  16  uniform  elements,  and  a  more  detailed 
one  (8x8)  containing  64  uniform  elements. 


Free  Vibration  Analysis 

Each  of  the  three  finite  element  models,  a  complex  model  (1254 
bar  elements)  and  two  equivalent  plate  models  (4x4  and  8x8  meshes), 
was  analyzed  to  determine  the  natural  frequencies  and  mode  shapes  of 
the  structure.  The  lower  natural  frequencies  and  the  associated  mode 
shapes  predicted  by  the  three  models  are  shown  in  Fig.  6.  The  mode 
shapes  predicted  by  all  the  models  are  very  similar.  Contour  plots 
of  several  of  the  elastic  modes  generated  for  the  structure 
contained  in  Appendix  A.  The  three  sets  of  frequencies  are  in  good 
agreement,  especially  at  low  frequencies.  The  relative  costs  to 
execute  each  vibration  analysis  are  given  in  Table  1.  Even  with  as 
few  as  eight  cell  repetitions  in  each  direction  the  cost  savings  using 
the  equivalent  plate  finite  element  model  are  substantial. 
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□  COMPLEX  MODEL  (1254  DOFS) 
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Figure  6.  Elastic  Modes  of  the  Space 
Lattice  Structure. 
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Table  1.  Relative  Vibration  Analysis  Costs 


Model 

Relative  Cost 

Complex 

1.0000 

4x4  elements 

0.0350 

8x8  elements 

0.1000 

Control  Problem 

The  control  problem  used  to  evaluate  the  capability  of  the 

equivalent  continuum  models  involves  angular  sensors  and  actuators 

located  at  the  center  and  the  four  corners  of  the  structure.  This  is 

similar  to  a  case  examined  for  an  isotropic  plate. ^  It  is  assumed 

that  the  output  is  measured  and  the  disturbance  is  acting  on  these  ten 

rotational  d.o.f.'s  at  the  five  locations.  Therefore,  the  [B],  [D], 

[C]T,  and  [M]T  matrices  are  identical.  The  output  and  control 

weighting  matrices  are  assumed  to  be  identity  matrices.  The  disturbance 

weighting  matrix  [W]  is  assumed  to  be  0.0001  times  the  identity  matrix. 

-12 

The  measurement  weight  matrix  [V]  is  assumed  to  be  10  times  the 
identity  matrix.  The  modal  damping  c  is  assumed  to  be  0.005. 

Evaluation  Model 

A  system  containing  24  modes  was  used  to  design  a  control  system 
using  LQG  theory.  The  evaluation  model  was  chosen  to  include  the 
two  rigid  body  modes  that  are  both  observable  and  controllable. 

The  frequencies  and  modal  coefficients  are  given  in  Table  2. 
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Table  2.  Evaluation  Model  Rotational  Coefficients 


Natural 
Freq 
(Hz)  . 

Rotation  about  the 

y-axis  (0X) 

Rotation  about  the 

x-axis  (8V) 

Corner 

1 

Corner 

2 

Corner 

3 

Comer 

4 

_ 

Center 

Corner 

1 

Corner 

2 

Corned 

Comer 

4 

Center 

i  - 

0.00 

-3.76 

-3.76 

-3.76 

-3.76 

-3.76 

-3.76 

-3.76 

-3.76 

-3.76 

-3.76 

0.00 

-3.55 

-3.55 

-3.55 

-3.55 

-3.55 

3.65 

3.55 

3.55 

3.55 

3.55 

1.16 

5.88 

5.88 

-5.88 

-5.88 

0.00 

5.88 

-5.88 

-5.88 

5.88 

0.00 

1.99 

6.98 

-6.98 

-7.04 

6.99 

0.00 

-6.98 

-6.98 

7.03 

6.98 

0.00 

2.26 

6.49 

-6.50 

-6.67 

6.53 

0.00 

6.51 

6.54 

-6.39 

-6.53 

0.00 

2.76 

-11.72 

9.35 

■11.97 

9.29 

1.45 

-7.57 

2.75 

-7.21 

2.75 

2.95 

2.76 

2.75 

-7.57 

2.75 

-7.57 

2.95 

-9.35 

11.71 

-9.35 

11.71 

-1.45 

4.14 

7.92 

-8.50 

-9.01 

7.67 

0.00 

8.19 

8.16 

-7.71 

-8.21 

0.00 

4.23 

10.75 

10.17 

10.79 

11.08 

-2.87 

.43 

.80 

-.74 

.16 

-.06 

4.23 

.47 

.12 

.47 

.12 

-2.87 

-10.61 

■10.60 

■10.61 

-10.60 

2.88 

4.67 

-9.54 

-9.63 

9.81 

9.49 

0.00 

8.82 

-8.68 

-9.40 

8.68 

0.00 

4.79 

-10.30  - 

10.39 

10.67 

10.29 

0.00 

-10.92 

11.06 

10.30 

•11.06 

0.00 

5.84 

12.75 

9.94 

13.47 

9.43 

3.27 

8.32 

-1.90 

7.04 

.59 

.94 

5.85 

2.02 

-8.29 

2.02 

-8.29 

-.90 

9.88 

12.74 

9.88 

12.74 

3.26 

6.30 

-6.47 

6.03 

7.00 

-7.02 

0.00 

5.81 

6.19 

-7.25 

-5.92 

0.00 

6.42 

-5.78 

5.58 

6.25 

-6.19 

0.00 

-6.44 

-6.19 

5.65 

6.39 

0.00 

6.76 

-11.80 

12.83 

■12.53 

12.38 

.72 

2.57 

-4.29 

4.05 

-4.22 

-2.61 

6.78 

-4.26 

2.29 

-4.26 

2.29 

0.00 

-12.51 

11.99 

■12.51 

11.99 

0.00 

7.31 

-10.13 

-9.36 

10.04 

9.75 

0.00 

-9.76 

9.56 

8.97 

-9.72 

0.00 

7.62 

-12.01 

12.32 

12.20 

■11.91 

0.00 

-10.46 

f-10.96 

9.30 

10.76 

0.00 

7.68 

4.03 

-4.35 

-4.33 

3.86 

-.05 

-7.23 

-6.68 

8.65 

6.90 

.08 

8.38 

6.67 

7.37 

7.21 

8.43 

4.20 

-.81 

-4.70 

-3.07 

-4.24 

-1.71 

8.40 

-4.51 

-1.13 

-4.51 

-1.13 

-1.64 

-7.93 

-6.62 

-7.93 

-6.62 

-4.23 

8.81 

4.81 

9.41 

-9.20 

-4.69 

.22 

-7.12 

3.80 

7.19 

-7.50 

.75 

Modal  Cost  Analysis 

The  modal  costs  are  calculated  for  each  of  the  three  analysis 
models  and  plotted  versus  frequency  in  Fig.  7.  The  8x8  element  model 
and  the  complex  model  have  very  similar  modal  cost  distributions. 

Not  surprisingly,  the  4x4  element  does  not  predict  the  modal  costs 
as  closely.  The  frequencies  and  modal  coefficients  for  the  two 
finite  element  modal  models  are  given  in  Tables  3  and  4  for  the 
4x4  and  8x8  meshes,  respectively. 
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Control  Design  Models 

The  reduced  models  chosen  contain  the  two  observable  and  control¬ 
lable  rigid  body  modes  plus  up  to  ten  elastic  modes.  The  elastic 
modes  retained  are  those  with  the  largest  modal  costs.  The  performance 
of  several  reduced  order  controllers  is  shown  in  Fig.  8. 

Concl usions 

A  basic  procedure  to  analyze  plate-like  space  lattice  structures 
using  finite  element  models  of  equivalent  continuum  formulations  of 
the  strain  energy  and  kinetic  energy  of  the  structure  has  been 
developed  and  implemeneted.  By  retaining  the  transverse  shear 
deformations  in  the  finite  element  formulation,  the  equivalent  plate 
finite  element  models  are  found  to  give  modal  results  consistent  with 
those  generated  from  a  detailed  truss  bar  element  model  of  the  space 
lattice  structure  even  for  small  numbers  of  unit  cell  repetitions. 

The  modal  models  resulting  from  the  finite  element  model  are 
shown  to  serve  as  accurately,  yet  are  more  simplified  in  the  control 
design  process  compared  to  the  complex  truss  bar  element  model  for 
the  space-lattice  structure.  This  study  has  shown  that  choosing 
control  design  models  based  on  modal  sequence  may  involve  unnecessary 
modes.  By  using  Modal  Cost  Analysis  to  consider  the  control  problem 
based  on  a  sequence  of  lower  modes,  it  is  found  that  several 
intermediate  modes  do  not  contribute  much  to  the  control  design. 


I?  MODF  EVAI  U  AT  ION  MODEL 
(ALL  CURVES) 

DYNAMIC  ANALYSIS  METHOD 

o COMPLEX  MODEL 
a8x8  ELEMENTS 
□  4x4  ELEMENTS 


Q 

©-  -a 


CONTROL  ENERGY  Vu  (nt2-m2) 


Figure  8.  System  Performance  of  the 
Example  Lattice  Structure. 


CHAPTER  III  -  EQUIVALENT  MODELLING 
AS  A  DESIGN  TOOL  -  PARAMETRIC  STUDY 

Lattice  plate  finite  elements  based  on  a  continuum  model  of  a 

large  plate-like  lattice  space  structure  are  used  to  examine  the 

effect  of  variation  of  several  fundamental  structural  parameters  on 

18 

the  natural  frequencies  and  mode  shapes  of  the  structure.  Reduced 
order  controller  design  models  are  developed  using  modal  cost  analysis 
to  rank  the  modes  for  each  set  of  structural  parameter  values.  The 
linear  quadratic  Gaussian  (LQG)  controller  design  method  is  used  to 
develop  feedback  control  systems  for  each  set  of  structural  parameter 
values.  The  resulting  system  performance  is  then  evaluated  by  examining 
the  steady  state  regulation  cost  of  the  structure  as  a  function  of  the 
structural  design  parameters. 

Problem  Statement 

The  initial  space  lattice  geometry  used  in  this  study  was  used 
in  Chapter  II  to  demonstrate  the  effectiveness  of  equivalent  finite 
element  models  in  designing  feedback  control  systems  for  plate-like 
lattice  space  structures.  A  simple  unit  cell  (Fig.  9)  is  repeated 
eight  times  in  two  orthogonal  directions  to  generate  the  lattice 
space  structure.  The  lattice  space  structure  geometry  is  defined  by 
six  parameters.  Four  of  these  are  the  cross  sectional  areas  of: 
the  upper  surface  members  (A^),  the  lower  surface  members  the 


ORIGINAL 

VALUES 

-  UPPER  SURFACE  MEMBERS  A,=80mm2 

-  LOWER  SURFACE  MEMBERS  A^Omm2 

-  VERTICAL  MEMBERS  A^IOmm2 

-  DIAGONAL  MEMBERS  Arf=  10mm2 

L=h=75m 


Figure  9.  Space  Lattice  Structure  Unit  Cell. 


vertical  members  connecting  the  two  surfaces  (Av),  and  the  diagonal 
members  connecting  the  two  surfaces  (Ad).  The  remaining  two  parameters 
are  the  thickness  or  height  (h)  and  half-width  (L)  of  the  unit  cell. 

For  simplicity,  it  is  assumed  in  this  study  that  all  members  are  of 
circular  hollow  cross  sections  with  outer  diameter  (dQ)  and  inner 
diameter  (d.).  The  ratio  between  the  two  diameters  (dQ/d^)  is  defined 
as  a.  The  plate-like  lattice  structure  is  assumed  to  be  free  along 
all  four  edges.  The  control  system  configuration,  shown  as  a  midplane 
in  Fig.  10,  consists  of  a  set  of  four  pairs  of  torque  actuators  (A) 
about  the  x  and  y  axes  respectively  located  along  the  diagonals  of 
the  plate-like  structure  at  a  distance  of  21.3  meters  from  each 
corner.  The  eight  actuators  are  driven  by  a  feedback  control  system 
using  eight  sensors  (M)  measuring  the  angular  motion  of  the  structure. 
These  angular  sensors  are  located  in  pairs  at  the  four  corners  of  the 
structure.  Five  pairs  of  disturbances  (D)  are  used:  four  pairs  are 
at  the  same  location  as  the  actuators  and  one  pair  is  at  a  distance 
of  21.3  meters  from  the  center  of  the  structure  along  a  diagonal. 

Cross  Sectional  Area  Variation 
The  first  parametric  variation  performed  in  this  study  is 
optimizing  the  cross  sectional  areas  (A^,  ^  Ay  and  Ad)  for  fixed 
lengths  (L  and  h).  Examining  the  strain  energy  terms  calculated  for 
the  equivalent  continuum  formulation  (see  Table  2  of  Ref.  7),  it  is 
apparent  that  the  bending  stiffnesses  depend  on  the  sum  of  A^  and  A2 
and  the  transverse  shear  stiffness  depends  on  Ad-  The  natural 
frequencies  and  mode  shapes  for  small  amplitude  transverse  vibration 
of  the  lattice  are  virtually  independent  of  the  in-plane  stiffness 


which  depends  on  the  difference  between  and  A2  (see  Table  2  of 
Ref.  7).  Therefore,  the  natural  frequencies  and  mode  shapes  for 
transverse  vibration  of  the  lattice  are  independent  of  the  relative 
magnitude  of  Ax  and  A2.  For  simplicity  of  this  study,  the  cross 
sectional  areas  (A$)  of  all  the  members  in  the  upper  and  lower 
surfaces  is  assumed  to  be  the  same. 

A$  =  Ax  =  A2  (III.l) 

The  cross  sectional  area  Ay  is  held  constant  at  its  original  value. 

The  total  mass  of  the  structure,  which  is  held  constant,  is  given  by: 

M  =  [(1  +  — )  (Ax+A2)  +|pv  +r  Ad]  (III. 2) 

where  p  is  the  mass  density  of  the  structural  members  and  d  =  *4i^+L2. 
Substituting  the  assumed  parameter  values  into  Eq.  (III. 2)  yields 
the  following  relation  between  A$  and  A^. 

A$  =  (69.167  -  4.167xl05Ad)  x  10"6m2  (III. 3) 

The  members  of  the  lattice  structure  are  extremely  slender. 

Before  examining  the  effect  of  varying  A$  and  Ad  on  the  natural 
frequencies,  mode  shapes,  and  performance,  their  effect  on  the  Euler 
buckling  load  capacity  of  each  member  of  the  lattice  is  examined  and 
shown  in  Fig.  11.  Since  A2  is  related  to  Ad  by  Eq.  (III. 3),  only 
Ad  is  used  in  presenting  results.  The  buckling  load  of  each  tubular 
member  is  nondimensional ized  by  dividing  its  value  by  the  buckling 
load  of  the  initial  lower  surface  members 


Ad  (mm2) 


Figure  11.  Allowable  Buckling  Load  vs.  Diagonal 
Member  Area. 
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P0  =  -T^  (III. 4) 

where 

E  =  71.7x10; 

»(a4-l)d? 

'o  = - b4 -  ;  and 

d0  =  ad,  . 

Figure  12  shows  that  as  Ad  is  varied  to  either  extreme,  Ad  or  the 

inversely  related  A$  will  become  very  small  causing  the  members  to 

have  unacceptably  low  values  of  the  Euler  buckling  load.  For 

reference,  the  Euler  buckling  loads  of  the  surface  and  diagonal 

members  in  the  initial  configuration  are  shown  in  Fig.  11  as  dashed 

lines.  Figure  12  shows  the  effect  of  the  cross  section  area  Ad 

(and  the  inversely  related  A$)  on  the  natural  frequencies  of  various 

modes.  It  is  noted  that  modes  4,5  and  modes  7,8  are  double  modes. 

It  is  interesting  to  note  that  all  of  the  present  natural  frequencies 

for  the  ten  lowest  modes  approach  maxima  when  Ad  is  between  30  and 
2 

50  mm  .  Apparently,  the  trends  of  these  curves  result  from  the 
compensating  effects  of  three  factors:  the  bending  stiffness,  the 
transverse  shear  stiffness,  and  the  rotatory  inertia.  When  Ad 

decreases  (or  A$  increases)  both  the  bending  stiffness  and  the 

rotatory  inertia  increase,  whereas  the  transverse  shear  stiffness 
decreases.  Figure  13  shows  the  effect  of  Ad  on  the  modal  cost  for 

each  mode.  As  the  Ad  value  decreases,  the  modal  cost  of  the  first 

mode  becomes  increasingly  dominant,  whereas  the  modal  cost  curves  for 


several  other  modes  cross  each  other.  Such  crossings  indicate  that 
different  modes  are  retained  at  different  Ad  values  in  the  reduced 
order  feedback  control  system  design  models. 

For  each  set  of  parameter  values,  a  series  of  reduced  order 
controllers  are  designed  using  the  LQG  method.  The  system  performance 
is  evaluated  using  the  full  evalation  model  driven  by  the  reduced 
order  controller.  A  typical  performance  curve  (Fig.  14)  relating 
regulation  cost  (V^)  to  control  energy  (V^)  is  generated  by  varying 
p  in  Eq.  (11.27).  The  performance  curves  for  all  of  the  cases  studied 
are  very  similar  and  are  contained  in  Appendix  B.  Figure  15  shows 
the  regulation  cost  (Vy)  versus  Ad  for  various  constant  values  of 
(solid  curves)  and  for  the  minimum  (circled  points).  These 
results  are  for  a  controller  designed  retaining  two  rigid  body 
modes  and  three  elastic  modes  in  the  reduced  order  controller  design 
model.  The  results  for  other  reduced  order  controllers  were  similar 

and  are  also  in  Appendix  B.  The  minimum  regulation  cost  is  obtained 

2 

for  Ad  values  between  10  and  30  mm  .  It  is  interesting  to  see  that 
designs  near  the  Ad  value  of  the  original  configuration  perform 
better  than  designs  near  the  peak  frequency. 

Thickness  or  Height  Variation 

As  a  second  parametric  study,  the  effect  of  varying  the  depth 
of  the  structure  (h)  is  examined.  The  variation  is  performed  such 
that  the  buckling  load  in  the  slender  truss  members  connecting  the 
surfaces  and  the  total  mass  of  the  structure  is  held  constant. 

For  tubular  members  with  a  constant  ratio  of  outer  and  inner 
diameters,  the  moment  of  inertia  (I)  is  proportional  to  the  cross 


sectional  area  (Ay).  Thus  the  Euler  buckling  load  in  the  vertical 
members  is  held  constant  by  maintaining  a  constant  ratio  between 
the  cross  sectional  area  (Av)  and  the  height  (h).  For  simplicity, 
the  cross  sectional  area  A^  is  assumed  equal  to  Ay. 


Av  "  Ad  "  (^xloA)  h 


(III. 5) 


Substituting  the  assumed  parameter  values  into  the  mass  Eq.  (III. 2) 
yields  the  relation  between  A$  and  h. 


/2 

where 

h‘  =  h/7.5m. 


As  h  increases,  the  nond i mens i ona 1 i zed  Euler  buckling  load  in  the 
surface  members  (Fig.  16)  decreases,  eventually  indicating  that 
these  members  are  too  slender.  The  natural  frequencies  of  the 
first  ten  elastic  modes  of  the  structure  are  shown  in  Fig.  17. 

Note  that  the  maximum  frequency  of  the  present  ten  modes  occurs 
at  thicknesses  near  22.5  meters.  The  trend  in  these  curves  is  caused 
by  the  compensating  effects  of  Ag  and  h  on  the  bending  stiffness, 
transverse  shear  stiffness,  and  rotatory  inertia.  As  h  increases, 

As  increases  Eq.  (III. 6).  Therefore,  the  increase  in  stiffness 
due  to  increasing  h  is  offset  by  the  reduction  in  stiffness  due  to 
decreasing  A  . 


The  modal  cost  as  a  function  of  thickness  (h)  for  the  various 

modes  is  shown  in  Fig.  18.  Since  the  modal  cost  curves  show  some 

crossings,  different  elastic  modes  will  be  retained  in  the  reduced 

order  models  for  different  ranges  of  thickness  (h).  Figure  19  shows 

the  regulation  cost  (V^)  versus  thickness  for  various  constant 

values  of  p  (solid  curves)  and  for  the  minimum  value  obtainable 

(circled  points).  The  reduced  order  controllers  for  thicknesses 
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less  than  15.0  meters  at  p  =  1.0x10  were  unstable.  The  performance 
results  are  for  LQG  controllers  based  on  two  rigid  body  modes  and 
three  elastic  modes.  Similar  results  were  obtained  for  reduced 
order  controllers  of  other  orders  (see  Appendix  B).  It  is 
interesting  to  point  out  that  for  this  example  the  optimal  performance 
always  occurs  when  the  thickness  (h)  is  near  7.5  meters,  whereas  a 
local  maxima  occurs  when  h  is  near  15  meters. 

A  procedure  has  been  demonstrated  to  utilize  equivalent  continuum 
finite  element  modelling  methods  to  efficiently  examine  the  effects 
of  parametric  variation  of  the  cross  sections  of  plate-like  space 
lattice  structures.  Using  a  specific  space  lattice  structural 
configuration,  two  sets  of  parameter  variations  were  performed.  In 
the  first  case,  the  length  of  all  the  members  as  well  as  the  mass 
of  the  structure  were  held  constant.  This  caused  the  area  of  the 
surface  members  to  be  a  function  of  the  area  of  the  internal  diagonal 
members.  The  second  case  involved  the  variation  in  the  thickness 
of  the  space  lattice  structure.  Again,  the  length  and  total  mass 
of  the  structure  were  held  constant.  The  allowable  buckling  load 
of  the  slender  internal  members  was  also  held  constant.  Thus  the 


cross  sectional  area  of  the  surface  and  internal  members  was 
determined  to  be  a  function  of  the  thickness.  For  both  cases, 
maxima  were  obtained  in  the  natural  frequencies  and  optimal 
performance  values  of  the  design  parameter  were  calculated. 
Modelling  the  structure  using  the  finite  element  continuum  model 
is  seen  to  allow  more  flexibility  in  treating  geometry,  boundary 
conditions,  attachments,  and  other  structural  complexities. 


CHAPTER  IV  -  EQUIVALENT  MODELLING  OF 
FRAME  STUCTURES  WITH  RIGID  JOINTS 


In  this  chapter  a  micropolar  plate  continuum  model  of  plate-like 

space  lattices  whose  frame  members  are  connected  with  rigid  joints 

is  derived.  The  resulting  continuum  model  is  used  to  derive  an 

equivalent  plate  finite  element  including  micropolar  rotations  and 

transverse  shear  deformations  as  nodal  d.o.f.'s  in  addition  to 
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displacements  and  displacement  derivatives. 

A  hexahedral  space  lattice  structure  with  rigidly  connected 
frame  members  is  used  as  a  demonstration  problem.  Several  sets  of 
frame  member  cross  sectional  properties  are  used.  The  exemplary 
plate-like  space  lattice  structure  is  assumed  as  rectangular  in  shape 
with  rigid  joints  and  with  four  repeating  cells  along  each  edge.  The 
natural  frequencies  and  mode  shapes  of  this  structure  with  all  edges 
free  are  calculated  with  and  without  the  micropolar  strains  and 
rotations  using  the  micropolar  plate  finite  element.  The  natural 
frequencies  and  mode  shapes  are  also  calculated  using  a  detailed 
finite  element  model  (with  every  structural  member  modelled  by  a 
beam-column  element)  of  the  plate-like  space  lattice  structure.  The 
natural  frequencies  and  mode  shapes  obtained  using  the  micropolar 
theory  are  shown  to  be  in  good  agreement  with  those  obtained  using 
the  detailed  finite  element  modelling,  for  the  present  example.  The 


natural  frequencies  and  mode  shapes  calculated  neglecting  the  micro- 
polar  terms  are  found  to  differ  substantially  from  those  obtained 
by  the  detailed  finite  element  model.  A  static  analysis  is  also 
performed  for  a  rectangular  plate-like  space  lattice  structure  with 
rigid  joints  and  with  eight  repeating  cells  along  each  edge.  The 
plate  structure  is  simply  supported  at  four  corners  and  loaded  at 
the  center  by  a  transverse  concentrated  load.  The  static  deflections 
of  the  structure  are  calculated  using  the  same  three  types  of 
modellings.  Again,  the  static  deflections  obtained  using  the  micro- 
polar  plate  finite  element  model  are  in  good  agreement  with  those 
from  the  detailed  finite  element  model,  while  neglecting  the  micropolar 
terms  causes  substantial  discrepancies. 

Equivalent  Continuum  Representation 
The  repetitive  nature  of  many  large  space  lattice  structures 
allows  them  to  be  modelled  as  an  equivalent  continuum.  The  smallest 
fundamental  building  block  of  such  a  structure  is  called  a  unit  cell. 
The  strain  energy  expression  of  the  unit  cell  is  obtained  using  a 
finite  element  model  with  each  frame  member  modelled  with  a  one 
dimensional  beam-column  element  containing  both  axial  and  flexural 
rigidity. 

ucell  =7(iI)T  rW"<} 

where 

Ucen  is  the  strain  energy  of  the  unit  cell; 


{q}'  =  ({qj}',  (q2)  ,  •••  {qn>‘); 

(q()T  -  («,,  vr  w,,  6Xi>  6Z.); 

(u-,  v..,  w. )  are  the  translations  in  the  (x,y,z)  directions 
at  node  i; 

(6  ,6  ,6  )  are  the  rotations  about  the  (x,y,z)  axes 

xi  yi  zi 

at  node  i; 

[Kcei i ^  1S  the  stiffness  matrix  of  the  unit  cell;  and 
n  is  the  number  of  joints  in  the  unit  cell. 

In  this  study,  the  stiffness  matrix  was  computed  using  the  NASA 

Structural  Analysis  (NASTRAN)  code,  with  each  frame  member  modelled 

14 

as  a  CBAR  type  element. 

The  derivation  which  follows  is  equally  valid  for  unit  cells 
containing  other  structural  elements  such  as  plates  and  membranes. 

The  translation  and  micro-rotations  in  the  unit  cell  are  assumed  to 
vary  linearly  through  the  thickness  (h)  of  the  plate-like  structure. 
The  orientations  of  the  various  translations  and  rotations  are 
expanded  as  a  first  order  Taylor's  series  in  the  z  direction  yielding 
ten  displacement  components  at  the  midplane  (see  Fig.  20). 

ui  =  u<xi-  *i>  *  zi  \<v  >i> 

vi  =  v(V  *i>  +  Z1  *y(xf,  V 


JOINT  TRANSLATIONS: 

MIDPLANE  TRANSLATIONS 


JOINT  ROTATIONS: 

z  MIDPLANE  MICROROTATIONS 


Figure  20.  Orientations  of  Translation,  Rotation, 
and  Microrotation. 
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W.  =  w{xis  y.)  +  zi  ez(xi,  y.) 

8yi  -  9x(*i’  *1>  +  zi  >*xz(V  *i>  <IV'2> 

6xf  ■  9y(xi-  >i>  +  2i  V(xr  yi) 

where  (u,v,w)  are  the  translations  in  the  (x,y,z)  directions; 

(<f>  ,  <t>)  are  the  rotations  about  the  (-y,  x)  axes;  (e  ,  e  )  are  the 
*  y  x  y 

microrotations  about  the  (-y,x)  axes;  £z  =  "gf  ;  pxz  =  dex/dz;  and 

MyZ  =  de^/dz  all  at  the  midplane  of  the  plate  continuum.  The 

rotation  (e  )  about  the  axis  normal  to  the  plate  is  ignored.  It  is 
zi 

noted  that  the  positive  sign  convention  of  each  rotation  is  based  on 
the  right  hand  rule.  Each  of  these  ten  displacement  components  at 
the  plate  midplane  is  expanded  as  a  second  order  Taylor's  series  in 
x  and  y  about  the  center  of  the  unit  cell.  For  example,  three  are 
shown  below. 

2  2 

u(x,y)  =  a2  +  xa2  +  ya3  +  ^-  a4  +  ^  a5  +  \  a6 

2  2 

v(x,y)  =  3y  +  X3g  +  yag  +  3  ajg  2  all  ^~2  a12  (IV. 3) 

2  2 

uyz(x,y)  -  a55  +  xa56  +  yag7  +  \  a58  +  ^  a59  +  \  a60 

Using  Eq.  (IV. 2)  and  Eq.  (IV. 3)  the  nodal  displacement  vector  { q}  can 
be  expressed  in  terms  of  the  displacement  Taylor's  series  coefficient 
vector  (a). 


{q}  =  [T]  {a> 
6nxl  60x1 


(IV. 4) 
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where 

CV  =  CCS]  Z.[S]] 


s.  0  0  0  0 

0  s.  0  0  0 

0  0  s.  0  0 

0  0  0  0  s. 

000  -s.  0 

0  0  0  0  0 

where 

^  X.2  x.y,  y.2 

[sf]  =  [l  «,  y,  -f  -f] 

Substituting  Eq.  (IV. 4)  into  Eq.  (IV. 1)  yields  the  strain  energy 
expression  of  the  unit  cell  in  terms  of  the  displacement  Taylor's 
series  coefficient  vector. 

ucen  =  1  ^a>TCT3TC<cell][T]{a}  (IV. 5b) 

The  fundamental  assumption  of  micropolar  continuum  theory  is 

that  the  local  rotation  at  a  point  (microrotation)  is  independent  of 
20 

the  translations.  The  strain  components  of  the  micropolar  continuum 


are  defined  in  terms  of  the  ten  midplane  displacement  components. 


classical  Dlate  strains: 


_  du  .  _  dv  .  0  _  du  .  dv 

exx  ix  ’  eyy  ~  dy  *  “xy  ~  ~5y  clx 

d<l>x  di  dc|>x 

Kxx  _  dx  ’  Kyy  ~  dy  *  ^Kxy  dy  + 


transverse  shear  strains: 


micropolar  strains: 


rx  =  1  "  ♦x*  '  V  ry  =  1  '  V  '  ey; 

.  d9X  .  _  .  _  1  ,d6X 

Jxx  dx  ’  pyy  "  dy  ’  yxy  “  1  '  dy  +  dx  ' 


These  strains  are  expanded  in  a  strain  Taylor's  series  in  x  and  y  about 
the  center  of  the  unit  cell. 


xz  xz  xz,x  xz.y'  xz,x  2  xz,xy  J  xz,y^  2 

2  2 
-0,o  „  .  o  w  .  o  x,  o  xy  ,  o  y 

yz  yz  yz,x  yz,/  yz.x^  2  yz,xy  2  yz,y^  2 


r  =  r° 
x  x 


r  =  r 

y  y 


(IV. 7) 


y  =  y°  +  y°  X  +  y°  y 
xx  ^xx  xx, x  xx ,y  J 


y  =  y°  +  y°  X  +  y°  y 

yy  yy  yy.x  yy.y  J 


_  0  ,  0  „  .  0 

y  =  y  +y  X  +  y  V 

xy  xy  xy,x  xy,y  J 


e  =  e°  +  c°  X  +  e°  V+  e°  9  ^  +  e°  ^  +  °  ,  Z 

z  fcz  z,x  z,r  z,x^  2  z,xy  2  z,y^  2 


y  =  y  +  y 
xz  xz 


2  ‘ 
0  v  .  °  . .  ,  o  x  ,  o  xy  ^  o  0  y‘ 

X  +  U  V  +  li  O  "TT~  +  U  -if-  +  u  / 

xz,x  xz,y  J  xz,x^  2  xz,xy  2  xz,y  2 


From  Eq.  (IV. 3),  EQ.  (IV. 6),  and  Eq.  (IV. 7)  we  can  generate  a 
relationship  between  the  strain  Taylor's  series  vector  {e}  and  the 
displacement  Taylor's  series  vector  {a). 


U)  =  [A  1  {a} 


(IV.8) 


,y2  =  a12;  en  =  exy,X  =  ?(u,xy}  +  \x2)  '  a 


(u.y2  +  V»xy^  "  a6  +  all;  e13  e  kxx,x  =  *x,x 


^  =  a  ■  c  =  ic®  =  a®  -  = 

<>xy  -  a35’  e15  -  Kvy,x  ^y.xy  "  a41* 


J>y2  =  a42;  c,7  =  K°y>x  =  ?(^>xy  +  *y>x2) 


a35+a40 


2  * 


t/,0  +  v  _  a36+a41  . 

-  x,y2  ^y.xy'  2  ’ 


:  a43;  c20  '  2exz  w,x  +  ^x  -  a14  +  a31 ’ 


\  +  -  ai^  +  a^7>  E??  =  2ev,  v  =  w°v2  +  * 


and 


r  =  o  _  l,o0  .  Oo  x  _  23  28  . 

e52  uxy,x  2'ex,xy  y  2  * 

_  Ao  .  fio  x  _  a24+a29 
e53  =  vxy,y°  '  6x,y2  +  6y,yxy)  "  2 

The  rigid  body  motions  of  the  structure  {b}  can  also  be  expressed 
in  terms  of  the  displacement  Taylor's  series  vector  {a}. 

(bl  =  [Ab]{al  ( 

where 

{b}^  =  (u°,  v°,  w °,  <J>°,  4>y >  0»  &)  » 

6  *7  <“°y  -  v!x);  and  2  (ex,y  *  ey,x>  ' 

Equation  (IV. 8)  is  merged  with  (IV. 9) 


I* 


Be40,17  =  2*0;  Be41,15  =  1,0;  Be42,16  =  1,0;  Be47,29  =  1,0; 


Be48,30  =  L0;  Be49 ,33  =  1,0;  Be50,34  =  1,0i  Be51,35  =  1,0; 


Je52,36  =  1,0;  Be53,37  =  1,0;  Be54,38  =  1'0;  Be55,39  =  1,0; 


Be56,40  =  1,0;  Be57,41  "  1,01  Be58,42  =  1,0i  Be59,43  =  1,0; 


Be60  44  =  anc*  ot^er  terms  of  B  and  Bfa  are  equal 


to  zero. 


Substituting  Eq.  (IV. 11)  into  Eq.  (IV. 5)  and  taking  advantage  of  the 
fact  that  the  strain  energy  is  independent  of  rigid  body  motion  of 
the  structures  gives  the  expression  for  the  strain  energy  of  a  unit 
cell  in  terms  of  the  strain  Taylor's  series  coefficient  vector  { e> . 


“cell  =  7  {=>T[K]{C} 


(IV. 12) 


where 


[Kc]  -  [BE]TtT]T[Kcen][T][BE] 


Strain  components  in  the  plane  connecting  two  unit  cells  must  have 
identical  values  to  maintain  compatabil ity.  This  is  satisfied  by 
taking  the  origin  of  the  Taylor's  series  expansions  at  the  center  of 
the  unit  cell  and  requiring  that  the  expansions  for  these  strain 
components  be  even,  insuring  that  these  strains  are  identical  on 
opposite  sides  of  the  cell.  Therefore,  the  odd  derivatives  of 


5W 


Using  Castigliano's  Theorem,  this  is  the  same  as  requiring  that  the 
row  of  [K  ]  corresponding  to  the  strain  gradient  times  {e}  be  equal 
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.•V 


to  zero.  [K  ]  is  partitioned  so  that  all  of  the  constrained  terms 
in  Eq.  (IV. 1)  fall  at  the  end  of  the  vector  {e}. 

.T  i — .. 


U 


cell 


K 

rr 

£r 

K 

ro 

eoe 

0 

0 

Kxx’ 

Kyy* 

(IV. 15) 


'xz ’  yz ■ 


xx  Myy 


u°  )  are  the  strain  terms  which  are  constrained;  and  {e  }T  =  (e°  y, 

Aj  U  AA y  A 

t  ,  ....  y  )  are  the  strain  terms  constrained  in  Eq.  (IV. 14). 
x  x  y  y  xy  9  y 

The  constraints,  Eq.  (IV. 14)  are  then  rewritten. 

{0}  =  [K  T  K  ]  W 


e  e 

r  o  oo 


e» 


(IV. 16) 


This  equation  is  rewritten  as  an  expression  for  { e0>  the  constrained 
strain  vector. 


=  -[K  ]_1  [K  ]  (e_ } 

o  e  e  r 

oo  ro 


(IV. 17) 


Substituting  Eq.  (IV. 17)  into  Eq.  (IV. 15)  yields  the  strain  energy 
expression  in  terms  of  the  retained  terms. 


"cell  =1'=/  CW  {£r> 


(IV. 18) 


where 


]  =  [K 


e 


]  -  [K  ][K  ]-*[K  ]T 

rr  ro  oo  ro 


13x13 


6 


It  is  assumed  that  [K  ]  divided  by  the  area  of  the  unit  cell  (Aceii) 
forms  the  constitutive  matrix  of  the  continuum  at  all  points 


U 


structure 


1 

1 


{er} 


[Keq] 

IT*-  Ur>  dA 
Acell  r 


(IV. 19) 


Area  of 
the  structure 


Equivalent  Finite  Element  Model 
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The  16  d.o.f.  rectangular  plate  element  is  modified  to  include 
in-plane,  transverse  shear,  and  micropolar  rotational  terms  yielding 
a  40  d.o.f.  element.  The  interpolation  functions  are  assumed  as 
linear  except  for  the  out-of-plane  displacement  which  is  assumed  as 
a  bicubic  Hermite  polynomial. 

u°  =  Cj  +  c2x  +  c3y  +  c^xy 


v°  =  c5  +  c6x  +  c?y  +  c8xy 


r\  o  9  *5 

w  =  c9  +  c10x  +  cny  +  c12x  +  c13xy  +  c14y  +  c15x 

+  clgx2y  +  c1?xy2  +  c18y3  +  cigx3y  +  c2Qx2y2  +  c21xy3 
+  c22x3y2  +  c23x2y3  +  c24x3y3  (IV. 20) 


c25  +  c26x  +  c27y  +  c28xy 


2eyz  =  c29  +  C30X  +  C31y  +  C32xy 


6x  =  c33  +  c34x  *  c35*  +  c36x^ 


e?  *  c37  +  c38x  +  ‘39*  +  c40xy 

Using  the  strain  definition,  Eqs.  (IV. 6)  and  the  interpolation 
functions  (IV. 20),  the  strain  vector  {er>  is  related  to  the 
interpolation  coefficient  vector  {c}. 

(er>  =  [C]{C> 

where 
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The  strain  energy  expression,  Eq.  (IV. 19)  for  the  structure  is  used 
to  derive  a  finite  element  representing  a  portion  of  the  structure. 
Substituting  Eq.  (18)  into  Eq.  (16)  and  limiting  the  integration  to 
a  specific  element  gives  the  strain  energy  expression  of  the  element. 


uel-7—“{c)T  [C]T[KpQ][C]dA{c}  (IV. 22) 

e  acell  eq 

Area  of 
the  element 


The  40  d.o.f.  of  the  element  can  be  expressed  in  terms  of  the 
interpolation  function  coefficients  { c}  by  using  Eq.  (17)  and  the 
x,y  location  of  each  node  point  of  the  element. 


and  i  refers  to  node  c  of  the  element. 

The  strain  energy  expression  of  the  element  is  then  expressed  in  terms 
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of  the  nodal  d.o.f . 's. 

(IV. 24) 

[C]T[Keq][C]dA[R]  . 

Area  of 
the  element 

The  same  lumped  mass  matrix  used  for  the  repetitive  space  lattice 
with  pinned  joints  Eq.  (11.25)  is  used  in  this  chapter. 

Results 

Space  Lattice  Structure  Configuration 

The  hexahedral  space  lattice  geometry  (Fig.  21)  is  used  as  a 

demonstration  problem.  The  space  lattice  is  generated  by  repeating 

a  unit  cell  (Fig.  22)  several  times  in  two  orthogonal  directions. 

The  initial  properties  of  the  structure  were  chosen  to  be  similar 

21  22 

to  several  box  truss  geometries  which  have  been  developed.  * 

The  length  of  each  main  structural  member,  which  is  equal  to  the 
thickness  of  the  space  lattice  structure,  is  4  meters.  The  main 
structural  members  are  6.1  mm  thick  Graphite  Epoxy  Tubes,  11.7  cm 
in  diameter.  The  diagonal  braces  are  aluminum  tubes,  5.657  meters 
long  with  a  thickness  of  0.5  mm  and  a  diameter  of  1.58  cm.  The 
natural  frequencies  and  mode  shapes  are  calculated  for  a  repetitive 
plate-like  space  lattice  structure  with  four  cells  in  each  direction. 
The  static  displacements  are  calculated  for  a  centrally  loaded 
plate-like  space  lattice  structure  with  eight  cells  in  each  direction. 


uel  =  Keitel* 
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Figure  22.  Hexahedral  Space  Lattice  Unit  Cell. 


Natural  Frequencies  and  Mode  Shapes 

The  natural  frequencies  and  mode  shapes  of  this  space  lattice, 

with  free  edge  conditions,  are  calculated  using  three  different  finite 

element  models.  The  most  accurate  solution  among  the  three  is  the 

"detailed"  finite  element  mode,  which  is  obtained  using  NASTRAN, 

with  each  frame  member  of  the  structure  modelled  as  a  single  CBAR 

type  finite  element.  Guyan  reduction  is  used  and  the  six  rigid  body 

modes,  are  calculated  using  the  NASTRAN  free  body  support  option, 

Eqs.  (II. 7)  to  (11.16).  The  flexible  natural  frequencies  and  mode 

shapes  are  shown  in  Fig.  23  as  the  solid  curves.  Four  sets  of 

frequencies  are  obtained  by  assuming  four  values  of  A/AQ,  where  the 

-3  2 

value  of  cross  sectional  area  AQ  =  2.235x10  m  for  the  main 
structural  members  is  used  as  a  basis  for  nondimensional ization. 

It  is  assumed  that  the  ratios  of  cross  sectional  areas  and  moments 
of  inertia  between  the  main  structural  members  and  the  diagonal 
braces  are  as  follows: 

A  I  t  d  A  2 

=  89.4  ;  and  —  =  —  —  (— )  . 

A  I  t  3  A„ 


Therefore  the  main  structural  member  properties  scale  as  follows: 


t  -  (£■) 


1/2 
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I  =  (ir) 


(IV. 25 


the  diagonal  brace  properties  also  scale  using  Eq.  (IV. 25).  The 
equivalent  micropol ar  plate  finite  element  developed  above  is  then 
used  to  calculate  the  natural  frequencies  and  mode  shapes  which  appear 
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as  the  dashed  curves  in  Fig.  23.  The  agreement  between  the  results 
obtained  by  the  micropolar  and  the  detailed  finite  element  models 
is  good,  especially  for  the  relative  small  number  (four)  of  cells 
assumed.  The  equivalent  micropolar  plate  finite  element  can  be 
reduced  to  a  transverse  shear  type  plate  finite  element  by  removing 
the  terms  of  [K  ]  associated  with  the  micropolar  strains  and  con¬ 
straining  ex  and  to  zero.  The  natural  frequencies  and  mode  shapes 
calculated  using  this  equivalent  plate  model  are  shown  as  the  dashed 
dotted  curve  in  Fig.  23.  It  should  be  noted  that  except  for  the 
first  flexible  mode  the  transverse  shear  continuum  plate  model 
frequencies  deviate  significantly  from  those  results  obtained  by 
the  micropolar  plate  finite  element  model  and  the  detailed  finite 
element  model.  These  discrepancies  are  small  for  very  flexible 
tube  members  (as  A/Aq  is  small)  and  they  increase  as  the  tube 
members  become  stiffer  (as  A/AQ  increases). 

Static  Analysis 

The  static  deflection  is  also  calculated  using  these  three 
different  finite  element  models.  The  detailed  model  described  above 
is  used  to  model  a  quadrant  of  a  statically  loaded  plate-like  space 
lattice  structure  with  a  central  concentrated  load  normal  to  the 
plane  of  the  structure.  The  structure  is  constrained  against 
translation  normal  to  its  plane  at  each  corner.  As  in  the  calculation 
of  the  natural  frequencies,  four  values  of  cross  section  areas 
(A/AQ)  are  used.  The  resulting  nondimensional  central  deflection 
(w/h)  is  shown  as  the  solid  curve  in  Fig.  24,  where  h  is  the  thickness 
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Figure  23.  Flexible  Natural  Frequencies 
and  Mode  Shapes. 
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Figure  24.  Nondimensional ized  Central  Deflections  of 
a  Corner  Supported,  Centrally  Loaded 
Space  Lattice. 


of  the  plate-like  lattice  structure  as  shown  in  Fig.  22.  The 
micropolar  finite  element  models,  with  and  without  the  micropolar 
rotational  terms,  are  also  used  to  calculate  the  static  deflections 
which  are  shown  in  Fig.  24  for  the  four  cross  sectional  areas  as  the 
dashed  and  dashed-dotted  curves,  respectively.  As  in  the  case  of  the 
natural  frequencies  the  results  obtained  using  the  micropolar 
finite  element  model  agree  well  with  those  obtained  using  the 
detailed  finite  element  model.  Both  sets  of  results  deviate 
significantly  from  those  obtained  by  the  finite  element  model  that 
omits  the  micropolar  rotational  degrees  of  freedom  as  the  structural 
members  become  stiffer  (as  A/AQ  increases). 

Conclusion 

A  method  has  been  developed  to  allow  repetitive  plate-like  space 
lattice  structures  with  rigid  joints  to  be  modelled  effectively  as 
an  equivalent  micropolar  plate  continuum.  A  micropolar  plate 
finite  element  has  been  developed  using  the  strain  energy  expression 
of  this  continuum  model.  This  allows  considerable  simplicity  and 
flexibility  in  tresting  geometry,  boundary  conditions,  attachments, 
and  other  structural  complexities. 

It  is  found  that  the  present  micropolar  plate  finite  elements 
can  accurately  predict  the  static  deflection  and  natural  frequencies 
of  the  plate-like  lattice  structures.  This  is  done  by  comparing  the 
results  with  those  obtained  using  a  detailed  finite  element  model 
where  each  structural  member  is  modelled  using  one  CBAR  (6  d.o.f.'s 
per  node  beam-column)  element  in  NASTRAN.  The  present  example  study 


shows  that  the  micropolar  plate  finite  element  is  considerably  more 
accurate  than  the  one  without  micropolar  rotation  but  with  transverse 
shear  d.o.f. 's. 

While  numerical  results  have  only  been  generated  for  structures 
built  with  frame  members  connected  by  rigid  joints,  the  method  would 
certainly  apply  to  a  general  class  of  repetitive  plate-like  structures 
with  either  rigid  or  pinned  joints.  The  present  unit  cell  modelling 
method  as  described  in  Eq.  (1)  could  certainly  be  generalized  to 
include  plate,  membrane,  or  other  types  of  finite  elements,  without 
changing  any  subsequent  portion  of  the  method.  Logical  next  steps 
would  include  incorporating  semirigid  joints  into  the  continuum 
model,  extension  of  the  method  to  model  member  and  joint  damping 
as  continuum  properties  for  forced  response  analysis,  including 
geometric  nonlinearities,  and  applying  this  finite  element  model 
for  feedback  control  system  design. 


CHAPTER  V 


CONCLUSIONS 


Further  Developments 

In  this  section  several  potential  extensions  of  this  work  will 
be  discussed.  In  general,  the  theoretical  basis  of  the  idea  is 
presented  without  supporting  numerical  examples. 

Semirigid  Joints 

The  first  extension  considered  is  the  application  of  the 
micropolar  plate  finite  elements  developed  in  Chapter  IV  to  model 
large  repetitive  plate-like  space  lattices  with  semirigid  joints.  In 
order  to  use  this  element  for  these  type  structures,  the  detailed 
finite  element  model  of  the  unit  cell  would  be  developed  with  each 
frame  member  modelled  as  a  single  beam-columm  type  element  as  before. 
However,  there  would  be  a  separate  nodal  point  at  each  joint  connected 
to  each  element.  The  value  of  the  translational  displacements  for 
each  of  these  nodes  would  be  equated. 
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uimi 
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vimi 


(V.l) 
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where 


( u . - »  v..,  w..)  are  the  translations  along  (x,y,z)  of 

I  J  1  J  I J 

the  jth  frame  member;  and 


m^  is  the  number  of  frame  members  of  the  unit  cell  which 


intersect  at  joint  i. 


Each  rotational  d.o.f.  would  be  connected  with  rotational  spring 

elements  of  zero  length  (Fig.  25).  The  stiffness  matrix  of  such  a 

spring  element  in  terms  of  the  torsional  stiffness  (J)  of  the  spring 
12 

is  given  by: 


r j 

[Krot.spr]  =  _0  j 


with  this  stiffness  matrix  connecting  the  rotational  d.o.f.  pairs 

(eY  ,  ex  ),  (e  ,  6U  ),  and  (e  ,  R  ),  where  (ey  ,  Bv  ,  B7  ) 
xi,  xi  •  yi  yi<  zi  zi,  xi,  yi,  zi.- 

J  J  ^  J  J  J  J 

are  the  rotational  d.o.f.  about  the  (x.y.z)  axes  at  joint  i  for  the 
jth  frame  and  j  varies  from  1  to  the  number  of  frame  members  of  the 
unit  cell  intersecting  at  joint  i  (m. ). 

The  equations  of  constraint  Eq.  (V.l)  would  be  used  to 
substantially  reduce  the  order  of  the  stiffness  matrix  [^cei i 3 
the  unit  cell.  First  the  d.o.f. 's  { q}  are  reorganized  into  three 
vectors  {qr>  containing  one  set  of  translation  and  rotation  d.o.f. 's 
at  each  joint,  {qQ}  containing  the  remaining  rotational  d.o.f. 's 
at  each  joint,  and  {qc>  conaining  the  remaining  translational 
d.o.f.'s  at  each  joint.  The  matrix  [Kcei-|]  is  similarly  reorganized 
and  partitioned. 
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The  equality  constraints  Eq.  (V.l)  can  be  rewritten 

{qc>  =  [G]{qr} 


(V.  3) 
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Substituting  constraint  Eq.  (V.4)  into  the  strain  energy  expression 

1  7 

Eq.  (V.5)  yields  an  equivalent  strain  energy  expression. 
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[Kro>  CW  +  le]T[Koc3T5  and  [R00]  =  [K00]  . 


For  the  micropolar  plate  finite  element  which  is  developed  in  Chapter 
IV  the  mass  associated  with  the  rotational  d.o.f.  at  the  joints  of 
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the  unit  cell  is  ignored  since  the  "lumped  mass"  description  of 
each  individual  frame  member  is  used.  This  element  is  shown  to 
yield  good  natural  frequencies  and  mode  shapes  using  a  specific 
space  lattice  geometry.  It  is  consistent  with  this  assumption  to 
utilize  static  condensation  Eq.  (1 1. 4)  to  reduce  the  strain  energy 
expression  to  contain  only  terms  associated  with  {qr }. 

ucell  4<VTr'W‘V  <v-7> 

where 

[Rcell]  =  t’W1  '  £|(ro^Roo^  ^ro^  * 

It  can  be  seen  from  Eq.  (II. 5)  that  the  Guyan  reduction  is  exact  if 
there  are  no  mass  terms  associated  with  the  omitted  d.o.f.'s 
([Mro]  =  0  and  [Mq0]  =  0).  Therefore,  the  reduction  in  Eq.  (V.7) 
causes  no  additional  error  beyond  that  generated  when  lumped  mass 
matrices  are  used  for  the  beam-column  elements.  It  should  be  noted 
that  the  strain  energy  expression  Eq.  (V.7)  is  equivalent  in  both 
order  and  form  to  Eq.  ( I V . 1 )  and  could  be  used  to  derive  a  micropolar 
plate  finite  element  for  repetitive  plate-like  space  lattices  with 
semirigid  joints.  These  results  would  be  compared  with  those 
obtained  with  a  detailed  finite  element  model  developed  for  the 
entire  structure  in  the  manner  described  above  for  the  unit  cell. 
Numerical  results  were  not  generated  using  this  method  since  the 
resulting  detailed  finite  element  model  is  larger  than  can  be 
effectively  analyzed  using  NASTRAN  as  installed  at  Purdue  University. 


Consistent  Mass  Matrices 

The  results  obtained  in  this  study  were  obtained  using  lumped 
mass  matrices  for  the  individual  frame  members.  Such  matrices  contain 
half  the  mass  of  the  member  at  the  diagonal  location  associated  with 
each  of  the  translation  d.o.f.  of  the  element.  This  ignores  any 
rotatory  inertia  effects  within  the  individual  frame  members.  The 
results  of  the  numerical  examples  in  Chapter  IV  show  that  the  lumped 
mass  assumption  for  the  individual  members  is  reasonable.  One 
obvious  extension  of  this  work  is  to  develop  the  mass  matrix  in 
terms  of  consistent  mass  matrices  for  each  frame  member.  One  way 
to  implement  this  would  be  to  assume  that  the  spatial  derivatives 
of  the  velocity  components  at  the  plate  midplane  are  zero,  implying 
that  the  velocity  components  are  only  functions  of  time  for  each 
unit  cell.  Using  this  assumption  the  kinetic  energy  expression  of 
the  unit  cell  in  terms  of  the  joint  velocities  and  angular  velocities 
can  be  transformed  to  a  kinetic  energy  expression  in  terms  of  the 
midplane  velocity  components.  This  has  been  done  for  the  beam-like 

g 

micropolar  continuum  representation.  This  kinetic  energy  expression 
could  be  integrated  over  each  element  to  yield  the  mass  matrix  of 
each  element. 

An  alternative  method  would  be  to  apply  the  consistent  mass 

matrix  method  (where  the  kinetic  energy  expression  is  used  to  generate 

the  mass  matrix  of  the  element  in  the  same  manner  that  the  strain 

energy  expression  is  used  to  generate  the  stiffness  matrix  of  the 

3 

element)  to  the  equivalent  model.  The  kinetic  energy  expression 
of  the  unit  cell  calculated  using  a  detailed  finite  element  model 


of  the  unit  cell  is  used  as  the  initial  kinetic  energy  expression 
for  the  unit  cel  1 . 

K-E'cell  '  7  1^>T  tMcel  <V-8) 

where  { q }  is  the  vector  of  nodal  velocities  of  the  unit  cell  and 
[Mc0i i 3  is  the  mass  matrix  of  the  unit  cell  from  the  detailed  finite 
element  model  of  the  unit  cell.  If  identical  transformations, 
reductions,  and  integrations  were  applied  to  this  expression  that 
were  used  in  Chapter  IV  with  the  strain-energy  expression  Eq.  (IV. 1), 
the  result  would  give  the  consistent  mass  matrix  of  each  micropolar 
plate  finite  element  [M^]  defined  as  an  analog  to  [Kg-]  ]  of  Eq.  (IV. 24) 

Joint  and  Member  Damping 

There  is  considerable  ongoing  research  to  develop  space  lattice 

frame  members  and  joints  with  significantly  larger  damping  than  that 

2T 

of  traditional  metal  components.  In  this  study  all  of  the  control 
system  results  have  been  generated  using  an  assumed  value  of  modal 
damping  (c  =  0.005).  For  structures  with  larger  damping  associated 
with  specific  motions  of  the  structure,  it  might  be  desirable  to 
include  damping  effects  in  the  continuum  model.  One  method  of 
doing  this  would  be  to  build  a  detailed  model  similar  to  that 
derived  in  the  paragraph  above  dealing  with  rigid  joints  except 
that  the  element  damping  matrices  would  be  assembled  and  condensed 
to  yield  a  system  damping  matrix  for  the  unit  cell.  This  unit  cell 
damping  matrix  could  be  substituted  for  the  unit  cell  stiffness 
matrix  in  Eq.  (IV. 1).  As  for  the  consistent  mass  matrix  in  the 
preceeding  paragraph,  the  same  procedure  used  to  transform,  reduce 


and  integrate  from  the  unit  cell  stiffness  matrix  to  the  micropolar 
plate  finite  element  stiffness  matrix,  could  be  used  to  calculate 
the  damping  matrix  of  each  micropolar  plate  finite  element.  Once 
assembled  this  would  yield  the  equation  of  motion  of  the  structure 
as; 

[ K] { q }  +  [ B] { q  >  +  [M]{q>  =  {P}  (V. 9) 

where 

[K]  is  the  stiffness  matrix  of  the  micropolar  plate 
finite  element  model; 

[B]  is  the  damping  matrix  of  the  micropolar  plate 
finite  element  model; 

[M]  is  the  mass  matrix  of  the  micropolar  plate 
finite  element;  and 

{q }  is  the  vector  of  nodal  displacement  d.o.f.'s  associated 
with  the  micropolar  plate  finite  element  model  of  the 
space  lattice  structure. 

This  equation  could  be  solved  to  yield  complex  eigenvalues  and 
eigenvectors,  which  could  be  used  in  a  similar  manner  to  that 
described  in  Chapter  II  for  feedback  control  system  design  and  system 
performance.  Alternatively,  structural  dynamic  response  could  be 
generated  using  direct  integration  analysis  methods. 


Application  to  Non-Lattice  Type  Structures 
While  the  equivalent  continuum  models  in  this  study  were 
specifically  developed  for  large  repetitive  space  lattice  structures, 
their  usefulness  should  extended  to  many  other  types  of  structures. 

As  an  example,  one  such  application  will  be  discussed.  Large 
mirrors  are  often  designed  as  a  surface  plate  which  can  be  coated, 
supported  by  a  backup  structure  which  is  often  made  up  of  web  type 
structure  (for  example,  see  Fig.  26).  Detailed  finite  element 
models  of  this  type  of  structure  can  be  too  large  and  cumbersome 
to  use  effectively.  As  an  alternative,  the  micropolar  plate  finite 
element  developed  in  Chapter  IV  could  be  calculated  from  the 
stiffness  matrix  of  the  unit  cell  based  on  a  detailed  model  of  the 
unit  cell  (such  as  in  Fig.  27)  and  used  to  solve  for  the  structural 
response  of  the  large  mirror. 

Summary 

A  basic  procedure  to  analyze  large  repetitive  plate-like 
space  lattice  structures  with  pinned  joints  using  finite  element 
models  of  equivalent  continuum  representations  of  the  strain  energy 
and  kinetic  energy  of  the  structure  has  been  developed  and  implemented 
By  retaining  the  transverse  shear  deformations  in  the  finite  element 
formulation,  the  equivalent  plate  finite  element  models  are  found 
to  give  modal  results  consistent  with  those  generated  from  a  detailed 
truss  bar  element  model  of  the  space  lattice  structure  even  for  small 
numbers  of  unit  cell  repetitions. 


The  modal  models  resulting  from  the  finite  element  model  are 
shown  to  serve  as  accurately,  yet  are  more  simplified  in  the  control 
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design  process  compared  to  the  complex  truss  bar  element  model  for 
a  specific  large  plate-like  repetitive  space  lattice  structure.  This 
study  has  shown  that  choosing  control  design  models  based  on  a  modal 
sequence  may  involve  unnecessary  modes.  By  using  Modal  Cost  Analysis 
to  consider  the  control  problem  based  on  a  sequence  of  lower  modes, 
it  is  found  that  several  intermediate  modes  do  not  contribute  much 
to  the  control  design. 

A  procedure  has  been  demonstrated  to  utilize  equivalent  continuum 
finite  element  modelling  methods  to  efficiently  examine  the  effects 
of  parametric  variation  of  the  cross  sections  of  large  repetitive 
plate-like  space  lattice  structures.  Using  a  specific  space  lattice 
structural  configuration,  two  sets  of  parameter  variations  were 
performed.  In  the  first  case,  the  length  of  all  the  members  as  well 
as  the  mass  of  the  structure  were  held  constant.  This  caused  the 
area  of  the  surface  members  to  be  a  function  of  the  area  of  the 
internal  diagonal  members.  The  second  case  involved  the  variation 
in  the  thickness  of  the  space  lattice  structure.  Again,  the  length 
and  total  mass  of  the  structure  were  held  constant.  The  allowable 
buckling  load  of  the  slender  internal  members  was  held  constant. 

Thus  the  cross  sectional  area  of  the  surface  and  internal  members 
was  determined  to  be  a  function  of  the  thickness.  For  both  cases, 
maxima  were  obtained  in  the  natural  frequencies  and  optimal 
performance  values  of  the  design  parameter  were  calculated.  Modelling 
the  structure  using  the  finite  element  continuum  model  is  seen  to 
allow  more  flexibility  in  treating  geometry,  boundary  conditions, 
attachments,  and  other  structural  complexities. 
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A  method  was  developed  to  allow  repetitive  plate-like  space 
lattice  structures  with  rigid  joints  to  be  modelled  effectively  as 
an  equivalent  micropolar  plate  continuum.  A  micropolar  plate  finite 
element  has  been  developed  using  the  strain  energy  expression  of 
this  continuum  model.  This  allows  considerable  simplicity  and  flexi¬ 
bility  in  treating  geometry,  boundary  conditions,  attachments,  and 
other  structural  complexities. 

It  was  found  that  the  present  micropolar  plate  finite  elements 
can  accurately  predict  the  static  deflection  and  natural  frequencies 
of  the  plate-like  lattice  structures.  This  was  done  by  comparing 
the  results  with  those  obtained  using  a  detailed  finite  element 
model  where  each  structural  member  was  modelled  using  one  CBAR 
(6  d.o.f.'s  per  node  beam-column)  element  in  NASTRAN.  The  present 
example  study  shows  that  the  micropolar  plate  finite  element  is 
considerably  more  accurate  than  the  one  without  micropolar  rotation 
but  with  transverse  shear  d.o.f.'s. 

Several  extensions  to  equivalent  plate  modelling  methods  have 
been  examined.  These  include  techniques  to  allow  the  existing 
equivalent  micropolar  plate  finite  element  modelling  method  to 
solve  structural  dynamics  problems  of  large  plate-like  repetitive 
space  lattices  with  semi  rigid  joints,  consistent  mass  matrices 
for  the  individual  frame  members,  and  element  damping  for  the  frame 
members  and  joints;  as  well  as  non-lattice  repetitive  plate-like 
structures  (such  as  mirrors).  The  equivalent  plate  modelling 
methods  could  also  be  extended  to  include  geometrically  nonlinear 
terms  allowing  buckling  of  the  structure  to  be  examined. 


Conclusion 


As  the  size  and  complexity  of  large  repetitive  space  lattice 
structure  increases  the  need  for  accurate,  efficient,  analysis 
methods  to  allow  the  structural  response  (both  static  and  dynamic) 
of  these  structures  to  be  predicted  efficiently  also  increases. 
Developing  detailed  finite  element  models  of  large,  complex 
structures  is  extremely  time  consuming  and  cumbersome.  Analyzing 
structural  response  for  such  complex  structures  is  also  extremely 
costly  and  for  sufficiently  large  and  complicated  structures  may 
be  impractical,  particularly  in  the  design  phase  where  many  analyses 
are  required.  Numerous  methods  have  been  developed  to  simplify 
and  automate  the  analysis  of  large,  complicated  structures.  While 
these  advances  are  crucial,  particularly  for  large  repetitive 
structures,  even  with  extremely  refined  software  to  manage  the 
generation,  assembly,  and  transfer  of  large  data  files  from 
different  portions  of  the  structure  to  some  system  model,  most 
of  these  methods  achieve  the  model  reduction  by  removing  internal 
d.o.f.'s  from  portions  of  the  structure  (substructures)  and 
representing  these  portions  of  the  structure  in  terms  of  the 
deformation  at  the  boundary  points  with  adjacent  substructures 
(retaining  sufficient  internal  d.o.f.  to  represent  the  loading  and 
internal  dynamics  of  the  substructure).  In  these  methods,  the 
resulting  model  of  the  system  is  tied  to  the  internal  geometry  of 
the  fundamental  structure.  In  general,  the  locations  where  node 
points  of  the  system  model  can  be  loaded  are  limited. 


As  an  alternative  to  potentially  cumbersome  detailed  modelling, 
numerous  studies  have  been  performed  where  a  standard  continuum 
finite  element  (such  as  a  plate  element)  is  used  to  model  a  discrete 
structure  by  using  an  equivalent  set  of  fundamental  properties 
(such  as  E,  v,  p,  t,  ...).  Several  studies  have  been  conducted  to 
develop  systematic  methods  for  determining  these  properties. 

Recently,  methods  have  been  developed  to  systematically  generate 
constitutive  matrices  of  equivalent  continuum  representations,  where 
higher  order  strain  terms  have  been  retained  consistent  with  the 
motion  the  unit  cell  of  the  repetitive  structure  is  capable  of. 

In  a  few  cases,  these  have  been  used  with  existing  finite  elements 
to  analyze  specific  space  lattice  geometries.  However,  in  general 
the  equations  of  motion  have  been  solved  in  closed  form  or  numerically, 
usually  for  simply  supported  structures.  To  the  authors  knowledge, 
no  systematic  effort  has  been  made  to  develop  higher  order  plate 
finite  elements  with  displacement  functions  specifically  chosen  to 
be  consistent  with  the  higher  order  strain  terms  in  these  equivalent 
continuum  representations. 

This  research  effort  has  extended  the  existing  equivalent 
continuum  representations  where  necessary  (no  such  representation 
existed  for  large  repetitive  space  lattices  with  rigid  joints). 
Rectangular  plate  type  finite  elements  have  been  developed  whose 
displacement  functions  are  consistent  with  the  equivalent  continuum 
representation  containing  transverse  shear  deformations,  with  and 
without  micropolar  strain  terms.  While  these  elements  are  intended 
to  model  structures  with  many  repeating  unit  cells  per  pi  ate- type 


finite  element,  it  is  gratifying  to  see  that  the  elements  accurately 
predict  the  response  of  the  current  examples  even  for  very  small 
numbers  of  unit  cells  per  element. 

The  equivalent  plate  finite  elements  representing  large 
repetitive  space  lattice  structures  have  many  advantages  compared 
to  more  conventional  methods.  One  significant  advantage  of  this 
method  over  detailed  modelling  using  a  substructuring  method  is 
that  the  number  of  and  location  of  nodal  points  required  to  model 
very  large,  complex,  repetitive  plate-like  space  lattices  is 
determined  by  the  number  and  accuracy  of  the  modes  required,  not  by 
the  internal  geometry  of  the  space  lattice.  Specifically,  there 
is  no  requirement  that  (or  advantage  to)  having  internal  node 
points  of  the  structure  located  at  unit  cell  boundaries  or  joint 
locations  in  the  lattice.  Large  repetitive  plate-like  space 
lattice  structures  of  the  type  considered  in  this  study  will  be 
used  as  space  platforms  to  which  will  be  attached  various  other 
structural  members.  While  solving  the  equations  of  motion  of  the 
equivalent  continuum  representation  of  a  free  floating  space 
platform  with  attachments  can  certainly  be  done  numerically,  it 
is  simpler  and  more  convenient  to  model  the  space  platform  as  an 
equivalent  higher  order  pi  ate- type  element  which  can  be  coupled 
with  other  structural  components  and  analyzed  in  the  free  floating 
configuration  using  existing  finite  element  solution  techniques. 

Some  space  lattices  will  be  generated  with  different  lattice 
configurations  in  different  regions  of  the  structure.  The  equiva¬ 
lent  plate  finite  element  formulation  applies  directly  to  this  type 


of  structure  by  simply  using  different  higher  order  constitutive 
matrices  in  different  regions  of  the  structure.  While  the  elements 
developed  in  this  study  were  rectangular  in  shape,  for  simplicity, 
elements  of  other  shape  could  be  developed  using  similar  methods. 
Thus,  plate-like  space  lattices  of  arbitrary  geometry  can  be 
readily  analyzed  using  this  method. 

This  research  represents  a  new  philosophy  to  the  structural 
analysis  of  discrete  structures.  That  is,  an  equivalent  continuum 
representation  of  suitable  order  to  account  for  the  possible  motion 
within  a  unit  cell  is  used  as  a  bridge  connecting  a  repetitive 
discrete  structure  with  a  consistent  order  discrete  finite  element 
model  of  the  structure.  This  method  gives  accurate  results  while 
greatly  simplifying  the  modelling  process  over  more  direct  methods 
of  going  from  the  discrete  structure  to  a  reduced  finite  element 
model  of  the  structure  for  large,  complex,  repetitive  structures. 
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APPENDIX  A 

Pin  Jointed  Truss  Mode  Shapes 
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Figure  A. 2.  Flexible  Mode  2  -  Contour  Plots 


APPENDIX  B 


Parametric  Study  Performance  Plots 
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Figure  8.1.  System  Performance  Plot  (k.  =  5  mm^). 
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Figure  B.2.  System  Performance  Plot  (A^  =  10  mm^). 
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Figure  B.4.  System  Performance  Plot  (Aj  =  50  mrn^). 
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Figure  B.6.  System  Performance  vs.  Diagonal 

Member  Area  (Controller  Area  =  12). 
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Figure  B.7.  System  Performance  Plot  (h  =  3.75m). 
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